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Abstract 

A detailed discussion of semiclassical trace formulae is presented and it is demonstrated how a 
regularized trace formula can be derived while dealing only with finite and convergent expressions. 
Furthermore, several applications of trace formula techniques to quantum chaos are reviewed. Then 
local spectral statistics, measuring correlations among finitely many eigenvalues, are reviewed and 
a detailed semiclassical analysis of the number variance is given. Thereafter the transition to global 
spectral statistics, taking correlations among infinitely many quantum energies into account, is 
discussed. It is emphasized that the resulting limit distributions depend on the way one passes to 
the global scale. A conjecture on the distribution of the fiuctuations of the spectral staircase is 
explained in this general context and evidence supporting the conjecture is discussed. 
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Introduction 

In a dynamical system deterministic chaos manifests itself in an effective unpredictability of the dy- 
namics for large times. Different degrees of a chaotic behaviour can be observed by identifying the 
properties of ergodicity, mixing, positive dynamical entropies, etc. [|l|. All of these concepts require 
to investigate the limit t — > oo of the dynamics. In the following we will restrict attention to the 
-important- cases of autonomous, bound Hamiltonian dynamics. Nevertheless, one can still find ex- 
amples for all types of different chaotic behaviour. However, upon quantizing such classical dynamical 
systems chaos in the above sense disappears. In a loose manner one could explain this observation 
by the regularizing effect the quantum mechanical uncertainty principle has on an irregular classical 
dynamics. More precisely, the discrete spectrum of the quantum Hamiltonian H enforces the time 
evolution in quantum mechanics to be almost periodic. In contrast, an almost periodic classical time 
evolution is only observed for integrable systems. If one prepares a quantum system at time to = in 
a state tpo, the solution -(/'(i) of the Schrodinger equation at time t > reads 

V'(t) = ^c„(/.„e-i^"* , (0.1) 

n 

where ipn and En denote the eigenvectors and eigenvalues of the quantum Hamiltonian, respectively: 
Hifn = Enifn- As f — > CO the state described by ( |0.1| ) fluctuates in a possibly very wild manner, but 
cannot approach zero. More importantly, the same observation holds for the correlation function 



which has the same structure as ( p.ll) . Thus no mixing can occur, nor can any of the other features 
conventionally characterising chaos be observed. However, in the classical limit, which is formally 
obtained by setting Planck's constant /i = 0, one recovers integrable as well as chaotic dynamics and 
all possible situations in between. The precise behaviour of the classical system is dictated by the 
Hamilton function, which generates the classical time evolution. The latter can produce correlation 
functions that range from quasi-periodic ones to exponentially decaying ones. The reason for this 
seemingly paradoxical finding is that the two limits involved, ?i — > and t — > cxd, do not commute. 

The field of quantum chaos |2|, y] now is concerned with a search for fingerprints the chaotic 
behaviour of its classical limit leaves on a quantum system. Since due to (0.1) the spectrum and the 



eigenfunctions of the quantum Hamiltonian completely determine the time evolution of a quantum 
system, it seems natural to investigate statistical properties of eigenvalues and eigenfunctions. The 
goal of the following lectures now is to explain in some detail an approach to eigenvalue statistics 
employing trace formulae. In the first part we introduce semiclassical trace formulae and discuss some 
of their applications, including a discussion of the use of zeta functions. The second part then is 
devoted to an application of trace formulae to eigenvalue statistics. A certain aspect of the latter, 
namely the approach to the global distribution of eigenvalues, will be the central theme of the third 
part. 

1 Introduction to Trace Formulae and Zeta Functions 

The basic idea behind semiclassical trace formulae is to take the trace of the quantum mechanical time 
evolution operator, thereby loosing all information on eigenfunctions, and to exploit the Fourier duality 
between the time dependent and the energy dependent picture. In order to deal with well behaved 
quantities, in the course of the required manipulations one has to employ a certain regularization 
procedure. Finally then, one can express certain sums over the spectrum of a quantum Hamiltonian by 
sums over periodic orbits of the corresponding classical dynamics. In certain cases one is furthermore 
able to determine explicitly the leading order contributions as ?i — > 0. 



To get an idea of the kind of quantum systems to which one can apply the following procedure we 
now list three of the more prominent types of examples: 

1. Schrodinger operators of the form H = —i^/S. + V{x), where here A denotes the Laplacian for 
M and V{x) is a suitable potential. If V{x) — > oo for |x| -^ oo sufficiently fast, it is ensured that 
H has a discrete spectrum. The corresponding classical dynamics on the phase space W^ x M*^ 
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is generated by the Hamilton function H{p,x) = 2^ + V{x). 

2. Quantum billiards: Let D C M"^ be a compact domain with a sufficiently well behaved boundary 
dD. Then H = —2;^ A, where now A denotes either the Dirichlet- or the Neumann-Laplacian 
for D. The corresponding classical dynamics on the phase space D x S"^^^ is that of a free motion 
inside D with elastic reflections from dD. 

3. Let M be a compact Riemannian manifold of dimension d. Then H = —2^ A, where now A is 
the Laplace-Beltrami operator for the Riemannian metric on M. The classical phase space is 
provided by the unit cotangent bundle over M, and the classical dynamics is that of the geodesic 
flow, i.e., the free motion along geodesies. 

When it comes to explicit formulae, we will in the following primarily adhere to case 1. However, 
possibly after some more or less obvious modifications, the results carry over to the other cases as 
well. 

In all of the above cases the spectrum of H is discrete and bounded from below. One can therefore 
add a suitable constant to H in order to render the quantum Hamiltonian non-negative. Thus 

< ^1 < ^2 < ^3 < ••• (1-1) 

provides the list of quantum energies in which we count each eigenvalue with its respective multiplicity. 

1.1 Semiclassical Trace Formulae 

The crudest semiclassical approximation to quantum mechanics can be summarized in the following 
rule: Each eigenstate of the quantum Hamiltonian corresponds to a cell of volume {2ttK)'^ in the classical 
phase space. Brought in a more mathematical form this statement yields the leading semiclassical 
asymptotics for the spectral staircase, 

NiE):=#{n; En<E}^^^l le{E-H{p,x))dpdx, h^O. (1.2) 

For the density of states this means 

d{E) = ^N{E)=Y,HE-En)^j^jJ6{E-H{p,x))dpdx, h^O. (1.3) 

The ultimate goal of a semiclassical trace formula now is to provide correction terms to the asymptotics 
on the r.h.s. of ( |1.2D and ( |1.3D . But since the spectral density d{E) is a highly singular object, with 
distributional singularities at the eigenvalues En, one has to employ some smoothing procedure. 
To start with let us consider the time evolution operator U{t) = e~ft *, so that 

iPit, x) = (U{t) Vo) (x) = J K{x, y- 1) My) dy (1.4) 

is the solution of the Schrodinger equation with initial condition ip{0,x) = 'il^o{x). The distributional 
kernel K{x, y; t) then obeys the Schrodinger equation 

d - \ 

ih— - H^ j K{x, y; t) = ih 6ix - y) 6{t) , (1.5) 
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with initial condition 



Ihn K{x,y;t) = 6{x -y) . (1.6) 



It can for t > be expanded in an orthonornial basis of eigenfunctions {ipn{x)} of H, 

K{x,y-t) = Y,^n{x)^Jy)e-i^-\ (1.7) 

n 

SO that the formal trace of the time evolution operator is given by Tr tJ{t) = ^ e~h "*. This object 
does in general not exist as a smooth function of t because it has singularities, the most obvious one 
being at i = 0. In order to regularize the trace consider a smooth function p{t) with a compact support 
and define the bounded operator 

/ + 00 
p{t)U{t)dt, (1.8) 

-oo 

which has a finite trace, since 

dxj dtp{t)K{x,x-t) = Y,j p{t)e-i'^-'dt = Y,p(^-^j ■ (1-9) 

Here p{E) denotes the Fourier transform of the function p{t). Since p{t) was required to be smooth and 
compactly supported, its Fourier transform is a Schwartz-class test function, i.e., p{E) and all of its 
derivatives decrease faster than any power. Once then the eigenvalues En behave as En ~ const, n", 
n — > oo, for some constant a, the sum on the r.h.s. of ( |1.9[ ) is finite. Therefore, the map p h-> Tr?7[p] 
is a tempered distribution. 

We now choose the test function p{t) eh^^ and hence observe 

Tr?7[p(t)ei^*] = f dx f "^ dt p{t) e^^' K{x,x;t) = J2p (^^r^) ■ (1-10) 

For convenience we rename the test functions in such a way that ^{E) := p{—E), and thus ip{t) = 
2-Kp{t). Then dTToD reads 

^ ^ {^^) = ^ / ^"^ f^ ^* ^^*^ ""''* ^^'^' '^' *^ • ^^-^^^ 

Here ^{E) and (^(i) are required to be smooth functions, and ip{t) shall have a compact support. The 
strategy now is to find a semiclassical approximation for the r.h.s., which then yields a semiclassical 
approximation for the spectral density d{E). To see this one notices that 

hd{E + he) = ^h5{E + he-En) = X]^(^~ ^"/T^ j ' ^^'^^^ 

n n ^ ^ 

SO that upon evaluating ( |1.12| ) on the test function ip{e) one obtains 

{hd{E + he), ^{e)) =Y,^ (^^^) ■ (1-13) 

Some of the general features of a semiclassical approximation for the r.h.s. of ( |1.11 ) can be revealed 
once one chooses an ansatz for the kernel of the time evolution operator. Let us therefore assume that 

K{x, y- 1) = j^ I e^'^^^**'^) an{x, y, t; q) dq , (1.14) 



where g G M is some auxiliary variable. In order to achieve the initial condition ( |1.6| ) for K{x, y; t) 
at t = one demands that 



{x,y,Ql]q) = {x-y)q and aa(x,y, 0; g) = 1 . 



(1.15) 



Inserting the ansatz ( 1.14 ) into the Schrodinger equation (L5) for t > then yields an equation for 
the phase (p and the amplitude an- Assuming that the amplitude allows for a formal expansion in 
powers of h, 

an{x, y, t; q) = V(i^)^ ak{x, y, t; q) , (1.16) 



fc>0 



one compares like powers of h. As a result, one obtains to lowest order in h the Hamilton-Jacobi 

equation ||^, |5| 

d(t){x,y,t;q) , j^ /^ , / , n n n (^ ^^\ 
— \-H{V^(f){x,y,t;q),x) =0 (1.17) 

for the phase. The higher orders yield transport equations for the k-th order coefficients a^. In 
principle, starting with the Hamilton-Jacobi equation ( p.. 17 ) one could now successively solve these 
equations order by order and thus construct the time evolution kernel ( [1.14 ) to any desired precision in 
h. Our ambition in what follows, however, will be restricted to obtain the leading order semiclassical 
asymptotics. 

A solution (j){x , y , t] q) of the Hamilton-Jacobi equation ( [1.17 ) with initial condition ( |1.15| ) turns 
out to be closely related to a generating function for a canonical transformation in phase space that 
describes the classical dynamics backwards in time. One introduces S{x, q, t) := (j){x, y, t; q) + yq, 
which obviously solves the Hamilton-Jacobi equation (1.17) with initial condition S{x,q,0) = xq. 
Then, if {q, y) is an initial condition for the solution of the classical equations of motion that reaches 
the phase space point (p, x) at time t with energy E, the function S{x, q, t) generates the canonical 
transformation (p, x) i-^ iq,y)- The formalism of canonical transformations and generating functions 
[Hjsl hence provides one with necessary conditions to be imposed on the phase appearing in the ansatz 



d 



V^(/)(x, y, t;q) =p and Vj^0(x, y, t; q) = -q 

{x,y,t;q) = -E and Vq(p{x,y,t;q) = , 



(1.18) 



whenever {q, y) is the initial condition for a classical trajectory with energy E in phase space that 
passes through (p, x) at time t. 

We are now in a position to insert the ansatz ( 1.14 ) into ( 1.11| ), 



E^ 



En-E 

h 



2^^27rhr I'^^l'^^l'^l ^(*) e^nM-'^''-^i)+m ^(,^)^- a,ix, x, t; q) . (1.19) 



Inspecting the r.h.s. of ( 1.19 ) one realizes that an application of the method of stationary phase |5|, g] 
to the multiple integral allows to obtain the asymptotic behaviour of this expression as ?i — > 0. To 
this end one has to identify the points {x,t,q) at which the total phase (j){x , x , t; q) + Et becomes 
stationary, i.e., where 



= [Vx(l){x,y,t;q) +Vy(l){x,y,t;q)] 

d 

= —(l)(x,x,t;q)+E, 

= Vq(j){x,x,t;q) . 



y=x 



(1.20) 



A comparison with (|l.lg| ) now immediately reveals that ( |l.20D picks out those triples {x,t,q) such 
that (q, x) is a point in phase space that is connected to itself by a solution of the classical equations 



of motion with energy E and initial condition (q, y) = (q, x) at to = 0, and {p, x) = {q, x) at t. Hence, 
exactly when t is a period of some classical periodic orbit with energy E, and {q, x) is a point in 
phase space on such an orbit, then {x,t,q) is a stationary point for the total phase in ( |1.19 ). On the 
contrary, if t does not correspond to any period of a periodic orbit, no stationary point in the integral 
over {x,q) occurs. If one therefore chooses the test function Lp{t) in such a way that it vanishes on 
all periods of classical periodic orbits with energy E, the method of stationary phase |5|, ^ yields the 
estimate 0{h'^), as h ^ in (|1.19D . 



Apart from terms of 0{h°°), hence all relevant contributions to (1-19) come from periodic orbits of 
the classical dynamics. In this context a distinguished role is played by the time t = 0, since then the 
whole hypersurface Qe '■= {ip,x); H{p,x) = E} of energy E in phase space is one huge manifold of 
stationary points. The leading order as ^ — > of its contribution to ( 1.19 ) can, however, be explicitly 
calculated. To this end let x(^) be a smooth function with x(t) = 1 in a small neighbourhood of t = 0, 
which vanishes outside a somewhat larger neighbourhood of t = 0. It shall in particular vanish on all 
periods of non-trivial classical periodic orbits with energy E. The trivial identity 1 = xi^) + [1 ~ x(^)] 



will then be introduced under the multiple integral in ( |1.19 ) such that the contribution 



2TT{27rhy 



f dx I dt /dgxW'^(i)e^t'^(^'^'*'«)+^*] Y,{ih)''ak{x,x,t;q) 



(1.21) 



of the stationary points (x, t = 0,q) is separated from the further stationary points with t > 0. If one 
now introduces polar coordinates for q, q = Xlo, X = \q\, \uj\ = 1, one can perform the integrals over 
the variables t and A by the method of stationary phase. As a result, ( |1.21 ) yields 



(27r/i) 



d-l 



27r 



{i + o{h)} , h^o 



(1.22) 



see [y, ch.l2] for a closely related problem. If one exploits the relation ( 1.13| ), the expression ( |l.22| ) 
indeed reproduces for the spectral density d{E) the leading asymptotic term as given in (|1.3| ). 

Although the leading asymptotic behaviour has already been determined, the real challenge is to 
calculate corrections to this. One therefore has to consider the remaining part 



1 

2^ 



dx / dt [l-x{t)](p{t)eh^^K{x,x;t) 



(1.23) 



of (|l.llD , which due to the above considerations is modulo terms of 0{h°°) completely dominated 
by the contributions of non-trivial classical periodic orbits. Since in ( 1.23| ) the test function (p{t) is 
required to be of compact support the time integral indeed extends only over a finite interval. It 
therefore suffices to know K{x, x; t) for a certain bounded range of t-values. This observation turns 
out to be essential for the further manipulations performed in ( 1.23| ) in order to obtain a semiclassical 
approximation. Namely, the ffist step in this direction consists of obtaining the leading semiclassical 
asymptotics for the time evolution kernel. In this context one naturally considers /i — > for fixed t. 



An application of the resulting asymptotic expression in ( 1.23 ), however, requires to use the latter on 
the whole range of integration. If it were not for the test function (p{t), this would involve t -^ oo. But 
the two limits /i — > with t fixed, and t — > cxd with h fixed, do not commute, as already mentioned in 
the Introduction. In mathematical terms this means that the semiclassical estimate is not uniform in 
t. This remark should stress the indispensible role played by the regularization procedure applied in 
the course of the present discussion. 

At least for small t, the calculation of the leading order semiclassical asymptotics for the time 
evolution kernel goes back to Pauli 0. It reads 



K{x,y;t) 



1 



{2mhy. 



E 



det 



d^R 



'"l^,y 



dxkdyi 



,jrR-,,^y(x,y;t)-i^u-,,^y {1 + 0{h)} , 



(1.24) 



where the sum extends over all solutions ^x,y of the classical equations of motion with boundary 
conditions 7a;,j/(0) = y and Jx,y{t) = x. Furthermore, the quantity R^^y{x,y;t) is defined as the 
integral of the Lagrangian L{x,x) = ^x^ — V{x) along 7a;, y, 



Ri:.,y{x^y^t) ■= I L{lx,yit'),'yx,y{t'))dt' 
Jo 



(1.25) 



Finally, z/^^ ^ denotes the number of points on the trajectory jx,y which are conjugate to the initial 
point y. For small t, v^^ ^ vanishes, and this is the case covered by Pauli. The modification required 
when passing through conjugate points as t increases was later provided by Gutzwiller B. 



Inserting the semiclassical approximation ( 1.24 ) for the time evolution kernel into ( 1.23 ) yields the 
expression 



2-K {2-iTih) 



^—^ fdxfdt [i-xit)]mYl 



det 



d^R 



7a;, H 



dxkdyi 



y=x 



(1.26) 



,j[R-tx,xi^,^;i)+Et]-i^u.,^.^ n _|_ 0(h)} 



Since in general an explicit evaluation of (1.26) is out of reach one employs the method of stationary 
phase to obtain the leading semiclassical asymptotics. In order to determine the stationary points of 
the phase R'y^^^{x, x; t) + Et one has to solve 



= V^.i?^^,,(x,y;t)|^^^+ Vyi?^,_^(x,y;t) 



\y=x 



From classical mechanics Q one recalls the identities 



d 



^xR^x,y {x, y;t) =p , VyR^^^y {x, y; t) = -q , -gl^x,y i^, y; t) 



E, 



'lx,y ; 



(1.27) 



(1.28) 



where q and p denote the momenta along ^x,y at the instants to = and t, respectively. E^^^^ then is 
the energy of this trajectory. Thus, ( |1.27 ) picks out those closed classical trajectories ^x,x that share 
identical initial and final momenta, q = p', hence these are periodic orbits with energy E. That way 
one recovers the previous observation, see ( |1.19| )-( [l.2[l| ), that all relevant contributions to (1.10) are 
due to classical periodic orbits of energy E. However, one now is in a position to calculate the leading 
semiclassical contribution to ( |1.11D explicitly. 

As a first step, consider the contribution of non-trivial periodic orbits to the regularized Green 
function, 

G^,{x,y-E):=—J [1 - x{t)] 0{t) K{x,y;t) ei""' dt . (1.29) 



In (1.26) one therefore disregards the integration over x and moreover reintroduces y ^ x. A well 
known calculation of the integral over t with the help of the method of stationary phase, see for 
example ||^, ^, then yields 



G^,{x,y;E) = J2Al,yix,y;E) ei^--(^'^'^)-*i'^-- {I + 0{h)} 



l^,y 



Here 



At.S^^y^E):- 



0(T. 



lx,y) 



2-K(2mhy 



Dyx,y{x,y;E) 



(1.30) 



(1.31) 



is an amplitude factor attached to each classical trajectory '-)x,y connecting y and x with fixed energy 
E. The time T^^ needed for this may depend on the particular trajectory. Moreover, 



D-y^,y{x,y;E) 



M7. 



det 


(- 


dxkdyi J 




Q2 


«7.. 



det 



9 Stx,y 

dxkdyi 



a^s. 



"Ix,y 



dyidE 



dxkdE 

O ^1x,y 

dE'^ 



V. 



1x,y 



Ix.y 
+ 1 






(1.32) 
(1.33) 



Finally, S^^y{x,y;E) := R-y^,y{x,y;T^^y) + ET^^^^ can also be expressed as the integral of pdx along 
the trajectory 'yx,y 

In a second step, the trace of Gpo will be calculated by integrating (|1.3[l|) over the diagonal x = y 
with respect to x. Once again the method of stationary phase will be applied, and to this end the 
condition of stationarity for the phase. 



i^xSy^^y {x, y; E) + VyS^^^^ (x, y; E) 



x=y 



p — q = 



(1.34) 



yields all points x on classical periodic orbits with energy E as solutions. These points are obviously 
not isolated, as it would be required by the simplest version of the method of stationary phase, since 
already a single periodic orbit 7 constitutes a one dimensional connected manifold of stationary points. 
However, in case the set of stationary points x devides into a series Mj, j = 1, 2, 3, . . ., of connected 
smooth manifolds of dimensions rrij < d in configuration space, one can introduce local coordinates 
in suitable neighbourhoods of the Mj's such that x >-^ {u,v), where u = (ui, . . . ,Um) parametrizes 
Mj . The method of stationary phase is then applied to the integral over the transversal coordinates 
V = {vi, . . . ,Vd-m)- As a result, one obtains that 



TrG^ 



j 



h- 



A 



M,e 



i^M, 



{l + 0{h)] 



(1.35) 



where Am- is an appropriate amplitude which is independent of /i, and Sm denotes the constant value 
of the action S{x,x;E) on Mj. If 7 is an isolated periodic orbit, the corresponding manifold Mj is 
one dimensional, so that it yields a contribution of 0{hP) to Tr Gpo- A /c-parameter family of periodic 



orbits leads to rrij = A; + 1 and hence contributes 0{h 2). In particular, in a classically integrable 



2 



system each invariant torus is of dimension d so that its corresponding contribution is 0{h' 

When the classical dynamics is such that all periodic orbits are isolated, and thus all connected 
manifolds Mj of stationary points are one dimensional, the amplitudes Am- can be calculated explicitly. 
The result being due to Gutzwiller can, e.g., be found in |^, ^ and reads 



IV GJ. 



2-K ^ 



I det (M^- 1)1^ 



{1 + Oih)} 



(1.36) 



Here the sum extends over all classical periodic orbits with energy E such that their periods T^ are 
contained in the support of (p{t). T^^ then denotes the primitive period corresponding to 7, i.e., the 
period of the primitive periodic orbit 7^ attached to 7. If 7 is a A:-fold traversal of 7p we write 7 = 7p 
and obviously find that T^ = kT^^. Notice that k can be both positive and negative, corresponding to 



traversals of the primitive orbit in both directions. The quantity Mj appearing in (1.36) denotes the 
monodromy matrix, or stability matrix, of the periodic orbit 7. It is given as a linearization of the 
Poincare recurrence map on a surface of section transversal to 7 in phase space. Collecting now the 
contributions (1.22) and ( 1.36| ) to ( |1.11 ) yields the regularized Gutzwiller Trace Formula (GTF) 

'E„,-E\ wo\{Qe) (^(0) 



E 



95 



h 



(27r/i) 



d-l 



27r 



{l + 0{h)] 



Vmo ^"^ |det(M^^-l)|2 

where the sum over all periodic orbits with energy E has been replaced by a sum over primitive 
orbits and their fc-fold repetitions. Exploiting the relation ( |1.13 ) now allows to obtain a periodic-orbit 
representation for the spectral density |^ |2|, 

m = ^ {1 + 0(^)1 + -nEt ^'" '°'P" " 'f:-^ i' + om) ^ (1.38) 

(2,fi)<i ""trti |det(A/4-l)|I 

At this place a remark seems to be in order. The test function ip(t) cuts off the sum over periodic 
orbits in ( 1.37| ) since it has a compact support, i.e., it vanishes outside a finite interval [Ti,r2]. In 



hyperbolic classical dynamical systems the number of periodic orbits with periods not exceeding T is 
finite and obeys the asymptotic law 



phtopi 



M{T) := #{7; < r^ < T} ~ — — , T ^ 00 . (1.39) 

Here htop > denotes the topological entropy of the classical dynamics on the energy shell Qe- 



Therefore the sum on the r.h.s. of ( 1.37 ) is actually of finite length. However, the corresponding sum 



in (1.38) is infinite, and due to the exponential proliferation (|1.39|) of the number of periodic orbits it 



indeed is divergent. But this divergence is to be expected because the sum approximates the spectral 
density which is a singular object. Thus, (|1.38|) should be understood as a formal relation whose 



actual meaning is provided by ( 1.37 ) 



A glance at ( |1.38| ) reveals that for each value E at which the spectral density shall be evaluated one 
has to determine the classical perodic orbits of energy E, their actions, stabilities, etc. In general this is 
a formidable task which renders an application of the trace formula ( 1.3§| ) for d{E) almost impossible. 



However, for a certain class of dynamical systems considerable simplifications emerge in that periodic 
orbits need only be calculated at a fixed reference energy Eq. A rather simple scaling relation then 
determines all required quantities at arbitrary values E of the energy. Such a mechanism to apply 
requires a mechanical similarity which allows to associate uniquely a periodic orbit at each value E, 
given one at Eq. Furthermore, the actions of periodic orbits have to be homogeneous functions of the 
energy, S-./{XE) = A"S'.y(S) for all A > 0. One can now fix Eq and obtain the energy dependence of 
S.y{E) as 

S^{E) = E^ E^^ SyiEo) . (1.40) 

It then proves useful to discuss the GTF in terms of the scaling variable E°'. Examples for systems with 
a mechanical similarity and homogeneous actions can be found among all the cases 1.-3. mentioned at 
the beginning of section |l[ Hamilton functions H{p, x) = ■^ + V{x) with scaling potentials, V{Xx) = 
X'^V{x) for all A > 0, yield homogeneous actions of degree " = ^ + ^- Billiards and geodesic flows on 
Riemannian manifolds always show the property ( 1.401) with a = ^^ since S^{E) = J pdx = \/2mElj, 



where ly denotes the geometric length of the periodic orbit 7. We remark that in quantum systems 
whose classical limit is scaling in the above sense, the semiclassical limit ^ — > obviously is completely 
equivalent to the limit E" — > 00. Once a > this in turn is equivalent to the high-energy limit 
E — > 00. An error term of the form 0{h ), H ^ 0, can therefore be replaced by 0{E~ ), E°' -^ 00. 

As an example, let us now discuss the GTF for quantum billiards with two degrees of freedom in 
some more detail. For simplicity, units will be chosen such that 2m = 1, and by abuse of notation we 
write p := +yE > so that S.y{E) = pl^. The spectral density shall then be expressed in terms of 
the momentum variable p, 

d{E) = Y,HE-E^) = Y,5[p^-pl) = ^^Y.^{p-p^)=:^^d{p). (1.41) 



The trace formula for the spectral density d{p) thus reads 



1 



d{p) = doip) + -,j2y: 



U^ cos i^kl 



np 



h'^np 



2 '^f^'yp ) 



Jp k=l 



det ( M^^ 



l + 0( - 
P 



CX) 



(1.42) 



where (io(p) denotes the analogue to the first term on the r.h.s. of (1.38) and provides a mean behaviour 
for the spectral density. If d = 2 the leading order asymptotics follows from ( 1.22 ). In addition, the 
subleading term for do{p) is also known, 



doip) 



A 



2TTh^ 



PT 



L 



+ 



P 



oo 



(1.43) 



where A denotes the area of the billiard domain, and L the length of its boundary. The negative and 
positive sign correspond to Dirichlet and Neumann boundary conditions for the Laplacian, respectively. 
The sum over periodic orbits in (|l.42| ) does not converge, as has already been discussed for the 
corresponding sum in (1.38). Again, ( |1.42| ) shall rather be viewed as a distributional relation that 
has to be evaluated on suitable test functions in order to yield an analogue of ( |1.37 ). Due to the 
scaling property on does, however, no longer need an external parameter fixing the energy at which 
the periodic orbit sum is to be evaluated. Let us therefore choose a smooth test function h{p), whose 
further properties will follow from a subsequent discussion. Since the variable p derives from the 
energy variable E = p^, h{p) should be chosen as an even function, h{p) = h{—p), so that it can also 
be defined for negative p. An evaluation of ( 1.43|) on such a test function thus leads to a regularized 
trace formula for billiards. 



n 



h{pn) 



POO -1 °° 

/ mp) h{p) dp + -Y,Y1 

•^0 ^, 1 — 1 



^y^Jh]C-¥ 



np -' 7. 



7p fc=l 



det 
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{i + o{h)} , n^o 



(1.44) 



where the following definition enters, 

F^[h]{u) : 
Since 



h{p) cos ipu /i^ j dp . 



cos I pu 



vr 



in case ji^^ is even this yields 



J" fc [h] (n) 



r/^7 



i-^k^X'^^ 



^7 

-1)^~ cos(pn) 

A7 ^1 
-1) 2 sin(pu) 



H-y even , 
fi^ odd , 
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2^ 



+00 



h{p) e*P" dp 



-i^kji- 



■ip 



9{u) 



(1.45) 



(1.46) 



(1.47) 



The version of the GTF that emerges in this situation hence reads 



/" ~ 1 

X] ^(^") = / dfiip) h{p) dp+- 



EE 

7p k=l 



CO 1 p 

— hp ^ 



^7p 
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h 



det 



(n 



{1 + 0(^)1 , ?i^0, (1.48) 



and was in this form given in [10|. 

The criteria that fix the class of test functions h{p) to be admitted in the trace formula ( |1.48| ) 
derive from the necessity that all terms entering ( |1.4^ ) be finite. First of all, the leading asymptotic 
behaviour ( 1.43 ) of the spectral density yields Weyl's law 



N{p) := # {n; < p„ < p} 

9 



A 



Airh? 



P 



00 



(1.49) 



upon integrating d{p) once. Since obviously N{pn) = n, ( 1.49[ ) can be inverted to observe pn ~ 
const.^/n, n — > oo. This requires the test function to obey h{p) = 0{\p\~'^~^), \p\ — > oo, for some 
5 > 0, in order to render the sum on the l.h.s. of ( |1.48| ) convergent. Due to ( 1.43| ) the same condition 
also ensures the convergence of the first term on the r.h.s. of (1.48). The exponential proliferation 
( |1.39| ) of the number of periodic orbits leads one to anticipate that the Fourier transform g{u) of the 
test function h{p) has to be required to decrease exponentially for \u\ — > oo, i.e., g{u) = ©(e"'*^^^-""') 
for some e > and some characteristic constant a > that is determined by the distribution of 
classical periodic orbits. Due to the definition (1.47) of g{u) as a Fourier transform the exponential 
asymptotic estimate for g{u) is clearly equivalent to demand that the test function h{p) itself be 
holomorphic in a strip |Imp| < a + e. 

A means to characterize the constant a that determines the strip of holomorphy to be demanded 
for the test function h{p) is provided by the thermodynamic formalism, see [^, |l^. In the latter 
theory one introduces the Ruelle zeta function 



Cp{s):='[\(l-e 



Sh'yj-, j.D*j,' 



■7p 



-1 



(1.50) 



Iv 



where e"^p > 1 denotes the modulus of one of the eigenvalues of the monodromy matrix M^^; the 
other eigenvalue of this matrix then is of modulus e~"^p < 1. Thus the product over primitive periodic 
orbits in ( |1.50| ) converges for those complex s such that Res is large enough; here /3 is considered as 
a parameter and is kept fixed. Indeed, one denotes the abscissa of convergence for ( LSOj) as P{fi), so 



that the condition Res > P{j3) defines the right half-plane where the Euler product in ( |1.5C ) converges 
absolutely. Thus, P{(3) has the following representation. 



P{I3) = inf <^ t G 



7p 



•■ip 



-fSu. 



''P < OO 



(1.51) 



In the thermodynamic formalism P(/3) is known as the topological pressure of the classical dynamics. 
Once the latter is hyperbolic, -P(/3) is shown to be of a certain universal form; this in particular implies 

l) 1^2 = e~2"Tp [1 + o(l)], the condition which ensures an 



that P(i) > 0. Now, since | det (M^^, 

:iod 



absolutely convergent sum over periodic orbits in (1.48) reads 



e 2 "fp 



7p 



»'f 



< oo 



(1.52) 

cxD, and thus one identifies 
contains only finite 



Hence, ( |1.52D is satisfied as long as g{u) = 0(e (^(2)+^)^!"!) for |n| 
a = hP{^)>0. 

In summary one concludes that the trace formula for billiards as given in ( 1 
quantities, once the test function h{p) satisfies the three conditions 

1. h{p) = h{-p), 

2. h{p) is holomorphic in the strip |Imp| < a + e for some e > 0, where a is a positive characteristic 
constant, a = hP{i2), 

3. h{p) = 0(|p|~^~ ) for \p\ — > CXD, where (5 > is arbitrary. 

At this point we would like to add a remark on the mathematical status of the trace formula. 
Gutzwiller originally derived [^, § his trace formula for the spectral density d{E), and obtained the 
relation ( 1.38| ). Slightly later, but seemingly completely independently, similar investigations were 
performed in the mathematical community. The first rigorous proof of a trace formula was given in 
a special case by Duistermaat and Guillemin [O]: The 'quantum Hamiltonian' they considered was 
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an elliptic positive pseudodifFerential operator of degree one on a compact smooth manifold without 
boundary; an example for such an operator is the square root of minus the Laplacian on the manifold. 
The corresponding classical dynamics then is generated by a 'Hamiltonian function' which is given 
by the square root of the kinetic energy for a single particle. Since thus the mechanical similarity as 
discussed above applies, one can choose Eq = 1 as a reference energy. Then the classical dynamics 
generated by the kinetic energy term and its square root, respectively, coincide on i}Eo=i- Translating 
now the result of [^] into a relation for the spectral density d{p), one recovers ( |1.42 ) with h = 1. 
Subsequently, generalizations in several directions were achieved. Finally, under certain assumptions 
on the class of quantum Hamiltonians and on the regularity of the classical dynamics, the semiclassical 
trace formula as given in ( 1.37] ) was proven [14, 15|. In principle, the strategy employed in |1^, 14, ^] 
was to start with an ansatz like ( [1.14 ) and then to give the integral on the r.h.s. as well as the expansion 
( |1.16| ) a precise mathematical meaning. Thereafter the multiple integral in ( 1.19| ) could be evaluated 
essentially by employing the method of stationary phase. In our presentation, we tried to discuss the 
trace formula in a manner that somehow interpolates between Gutzwiller's original investigation and 
the rigorous mathematical treatment, which itself requires an extensive technical apparatus. 



1.2 Some Applications of Trace Formulae 

In order to illustrate the use of semiclassical trace formulae in quantum chaos, we are now going to 
discuss some of their applications. For simplicity the following considerations will be restricted to the 
cases 2.-3. of the list at the beginning of section ffl. Hence the relevant form of the trace formula is 
given by ( 1.48 ). To begin with let us study the test function h{p) := e~^ *, t > 0, which yields on the 
l.h.s. of the trace formula the trace of the heat kernel, 



e^(t) = TYe-^* = j;e-^ 



-pit 



(1.53) 



For quantum billiards as well as for Laplacians on Riemannian manifolds the asymptotic behaviour of 
&fj{t) as t ^ 0, or equivalently as h^ 0, is well known to yield [l^ 



A 



N-l 



71=1 

with some appropriate coefficients Cn that can in principle be successively determined. The r.h.s. of 
the trace formula can be e^ 
is inserted in ( |1.48| ). Thus 



the trace formula can be evaluated once the Fourier transform g{u) = ,. e « of the test function 



/4^ 



'7p 



^ii^t) 



Mp) e-f'' dp + 






V^TrhH t^ ^ I det ( M^ - 1 ) I ^ 



7p 



7p 



(1.55) 



For /i — > each term in the sum over periodic orbits is of 0{h°°) so that in (1.54) all power-like 

contributions in h have to derive from do{p). A term- by-term inversion then yields the semiclassical 

expansion 

Ap 



do{p) 



L 



+ 0{h) . 



(1.56) 



27r^2 I ^^f^ 

As remarked earlier, the asymptotics for /i — > can be converted into an asymptotics for p -^ oo, so 
that ( 1.56 ) exactly reproduces ( 1.43 ). In fact, the above reasoning is a proper justification for (|1.43|) . 
Among the first applications of trace formulae in quantum chaos one finds efforts to introduce 
semiclassical quantization rules for classically chaotic systems. These should serve as substitutes for 
the semiclassical EBK-quantization scheme which only applies to classically integrable systems. A 
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first guess of a quantization rule could be to use the representation ( [1.42 ) for the spectral density d{p): 
one has to find classical periodic orbits and to calculate their lenghts, stabilities, and Maslov phases. 
Upon evaluating the r.h.s. of ( |1.42 ) with these data one would then approximate d{p) and thus be 
able to identify the singularities of the spectral density. However, due to the lack of convergence of 
the periodic orbit sum in ( |1.42| ) one has no good control over the quality of the approximation that 
occurs by cutting off the sum after finitely many terms. Therefore a regularization is called for which 
renders the periodic orbit sum absolutely convergent. One hence is advised to consult the smoothed 
trace formula ( 1.48| ) with a suitable test function. Having in mind to approximate the spectral density, 
as a natural choice for an admissible test function a Gaussian 



Kv') 



e^/^T 



^ + e ^ 



(1.57) 



was suggested in [17j, so that the l.h.s. of ( 1.4S| ) yields a smoothed spectral density. This approaches 
d{p) in the limit e — > 0, 



hni V h{pn) = V [6{pn -p)+ KPti + p)] = d{p) 



(1.58) 



if p > 0. The second Gaussian in ( |l.57| ) is necessary to yield an even test function, h{—p') = h[p'). 



If p is not too small, the first Gaussian yields peaks of height —^ at each p 



and the second 



Gaussian only adds a negligable value. Thus, upon scanning p one can detect the quantum energies as 

En = Pn- Now the r.h.s. of (1.48) can be evaluated with the Fourier transform g{u) = - cos{pu) e ^. 
Therefore the periodic orbit sum 
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EE 
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-i^kfi: 
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(1.59) 



clearly converges absolutely due to the Gaussian suppression of the exponential proliferation ( 1.39| ) 
of the number of terms entering ( 1.59| ). Moreover, one can explicitly observe that the limit e — > 
applied to ( 1.59 ) recovers the periodic orbit sum in ( [1.421) . 

Integrating the spectral density once with respect to p yields the spectral staircase ( [1.49| ). The 
trace formula with the test function ( 1.57[ ) will therefore now be integrated in p over the interval [0, q]. 
On the l.h.s. one then obtains the approximation 



0<pn<g 



E 



Pn=q 



- + 0[eer^ 



+ Eo 

Pn>q 



ee e^ 



(1.60) 



to the spectral staircase. However, in the limit e — > the expression ( 1.60| ) approaches the symmetrized 
staircase function 

Noiq) = lim i [N{q + 5) + N{q - 6)] , (1.61) 

which is identical to N{q) whenever q ^ pn- Otherwise one obvioulsy observes that N{pn) = n — i^. 
This allows to set up the quantization rule 



cos[7riVo(Q')] = <^ g e {pi,P2,P3,---} • 
According to (1.591) the leading periodic orbit contribution to ([1.6 Op reads 
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(1.62) 



(1.63) 
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which again is an absolutely convergent sum as long as e > 0. In 1 18] it was suggested to use ( 1.63 ) in 
( 1.62| ) to obtain an approximate semiclassical quantization rule. It was furthermore demonstrated in 
several examples that this procedure yields good numerical approximations to the quantum energies 
En = Pn^ and that this scheme is quantitatively superior to the direct use of a Gaussian test function 
( |1.57| ) in the trace formula. 

Inspecting the periodic orbit sum ( |1.63 ) one observes that, viewed as a function of q, it constitutes 



a superposition of sine-oscillations with wave lengths jj— . Adding those terms finally yields Ni;{q) 
One can now ask how many of these terms become efficient when one is interested in resolving all 
quantum energies En < p^ with the quantization rule ( |1.62|) , see |^] for related discussions. At the 
point p the mean density of zeros of the expression cos[i\'Nq{p)\ is given by ( 1.43| ) so that their mean 
separation Ap asymptotically reads Ap ~ ^5^- ^^ order to effectively resolve the zeros pn in the 



vicinity oi p with the help of the periodic orbit sum ( 1.63| ) one hence has to add all oscillatory terms 
with wave lengths down to Ap, i.e., one needs all periodic orbits 7 = 7p with |p^ > Ap - ^ 



7p " ^P 

Thus, a cut-off at approximately kl^^ ~ -^ ensures that one has taken all relevant contributions into 
account. Since the exponential proliferation ( |1.39| ) of the number of periodic orbits reads in terms of 
their geometrical lengths 



M{l) = *{r, l; 



<t}.- 



I 



OO 



(1.64) 



where r denotes a scaled topological entropy which is independent of p, the number Np of terms to 
be included in ( 1.63| ) in order to resolve all quantum energies with En < p^ is given by 
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Recalling now the asymptotics ( 1.4S| ), the number Nn of periodic orbits required to resolve the lowest 
n quantum energies thus derives from ( 1.65 ) to behave as 



Nn 
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2T\^TrAn 
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OO 



(1.66) 



Therefore, the computational effort to obtain more and more quantum energies by using the semiclas- 
sical quantization rule grows enormously. 

Similar considerations apply to all quantization procedures that are based on semiclassical trace 
formulae since the structure of the periodic orbit sums involved together with the exponential prolif- 
eration of the number of periodic orbits will always produce estimates of the effort to be spent which 
are similar to ( [L.66 ); one will certainly not achieve an improvement that comes down to a power-law 
behaviour for Nn- However, in some specific systems, especially in classically chaotic scattering sys- 
tems, it might happen that one can resolve a certain number of zeros using only a few periodic orbits. 
In these situations the specific features of the systems cause the number of required periodic orbits 
to blow up only at considerably large Pn- In this context one should bear in mind that ( 1.66| ) is only 



an asymptotic statement. The principal effect, resulting in the estimate ( 1.66 ), is however inherent 
in semiclassical quantization procedures based on trace formulae and cannot be simply overcome. At 
this point also the analogy between the zeros of the Riemann zeta function, see below, and quantum 
energies of classically chaotic systems breaks down since the density of Riemann zeros is considerably 
lower than ( |1.43p . This results in a power-law estimate for Nn-, and that has to be compared with 
( |1.66| ), see [[L9|] for a brief account. Basically the same problem occurred in the mathematical literature 
on the Selberg trace formula, where it was noted that well known estimates that hold in the theory 



of the Riemann zeta function could not be carried over to the Selberg zeta function, see [20| for an 
extensive discussion. 

As a consequence, various studies have revealed that it is possible to calculate a certain number, say 
some ten, quantum energies from classical data employing semiclassical trace formulae. However, the 
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huge numerical effort to be spent in order to increase the output of the number of eigenvalues prohibits 
to apply these semiclassical methods on a large scale. It turns out that numerical methods to solve 
the eigenvalue problem of the quantum Hamiltonian directly are far more efficient than semiclassical 
ones, if one is interested in obtaining a large number of quantum energies. 

1.3 The Selberg Trace Formula and the Selberg Zeta Function 

The trace formulae discussed so far essentially yield the leading semiclassical approximation to the 
spectral density. In addition, one has to cope with a further approximation in applications because in 
practise one can only take finitely many terms of a periodic orbit sum into account. Thus, it is often 
difficult to trace back numerical inaccuracies to either one of the approximations applied. Luckily 
there exist cases in which an exact trace formula is at hand, which therefore does not include any 
semiclassical approximation. The classical systems to be considered are geodesic motions on manifolds 
of constant negative curvature, and their quantizations are provided by the quantum Hamiltonians 

= — ^A, where A denotes the Laplace-Beltrami operator for the manifold. Due to the negative 
curvature the classical dynamics show a strong chaotic behaviour, and thus these models serve as 
convenient playgrounds for quantum chaos. 

In the simplest case of d = 2 degrees of freedom one has to deal with the geometry of hyperbolic 



surfaces, see for example [21, 22 1. Here we choose as a model for two dimensional hyperbolic geometry 
the unit disc 

V:={z = x + iy eC; \z\ < 1} (1.67) 

endowed with the Poincare metric ds'^ =4(1 — x^ — y^)^^ (dx^ + dy'^) of constant negative Gaussian 
curvature K = —1. The classical dynamics of a single particle on T> is that of a free motion along 
geodesies. The latter are those circular arcs and straight lines in D that are perpendicular to the unit 
circle. The hyperbolic distance d(z, w) between two points z ^ w inD is measured with the Poincare 
metric ds'^ and is given as the hyperbolic length of the unique geodesic arc connecting z and w. In 
explicit terms it reads 

2|2 — wl'^ 

cosh d{z,w) = 1 + — 1 .„, ,, ■. prr . (1.68) 

(1 — l^;] j(l — |if| j 

Prom this relation one immediately concludes that any point w on the unit circle is infinitely far away 
from any point z in the interior of D. 

Due to the latter observation the Poincare unit disc itself is an infinitely extended two dimensional 
hyperbolic space. In order to yield compact and closed surfaces, and thus examples for the case 3. 
of the list at the beginning of section [I], one has to confine oneself to bounded domains in the unit 
disc supplied with appropriate boundary identifications. In as much as a fiat torus can be constructed 
from a parallelogram in the euclidean plane with opposite edges identified, one can obtain compact 
surfaces with a metric of constant negative curvature by an analogous construction. The boundary 
identifications that produce a torus out of a parallelogram are translations on the euclidean plane. 
The group generated by the two identifications of a torus then is a discrete subgroup of the group of 
motions of the plane, i.e., the group of translations and rotations. The analogous transformations on 
the Poincare disc are the Mobius transformations 

az + p 
z ^ _ , (1.69) 

pz + a 

with |ap — |/9p = 1, which indeed map the interior of V to itself. Introducing the matrix group 
SU{1,1) := < ( ^— ) ; jap — |/3p = 1>, a composition of two Mobius transformations (1.69) corre- 



sponds to a matrix multiplication in SU{1, 1). The latter group is therefore the analogue of the group 
of motions on the euclidean plane. Given now a suitable bounded domain T <ZT), with finitely many 



geodesic arcs as its boundary components, the Mobius transformations ( 1.69 ) identifying pairs of edges 
generate a discrete subgroup P C SU{1, 1). 
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The quantum Hamiltonian H = —-^A requires to know the Laplace-Beltrami operator A for the 
Poincare disc P, 

If this operator is defined on the smooth functions on V it has a continuous spectrum [^,00). Due to the 
pecuhar value of the bottom of the spectrum one parametrizes the quantum energies as E = ^{p'^+j), 
with p > 0. For the fohowing we choose units such that h = 1 = 2m, and hence E = p"^ + ^. If one 
now constructs a compact and closed surface from a suitable domain T <Z T) as described above, 
one imposes periodic boundary conditions on the eigenfunctions of H = —A. Let gi,. . . ,gk be the 
Mobius transformations ( 1.69| ) identifying pairs of edges of J-, then one demands that il){gjz) = 4'{z), 



j = l,...,k. Since the transformations gi,---,gk generate the discrete group T C SU (1,1), the 
periodicity extends to all elements of T, ip{gz) = V'(-z) for all g GT. The Laplacian now being defined 



on a compact surface has a discrete spectrum = Eq < Ei < E2 < . ■ ■■ If one defines pn '■= \/ En — \, 
the spectral staircase satisfies the Weyl asymptotics (|l.49| ). 



In the 1950's Selberg |23| employed the above constructions to obtain the Selberg Trace Formula 
(STF) 

Y.hipn) = ^ rptanHnp)hip)dp + Y.fl ^""'^V ^ (1-71) 

n 2vr7o Vtt2sinh(^) 

see also [pO|. Here A denotes the area of the given surface, measured with the Poincare metric, and 
the sum over the 7p's extends over all primitive periodic orbits of the associated classical dynamics. 
Since the latter are given by the geodesic motion, the periodic orbits are the closed geodesies on the 
compact surface. The test function h{p) is required to fulfill 

1. h{-p) = h{p), 

2. h{p) is holomorphic in the strip |Imp| < 2 + £ for some e > 0, 

3. h{p) = 0{\p\^'^~ ) for \p\ — > 00, where (5 > is arbitrary. 

Furthermore, g{u) = -^ f^^ ^{p) e*^" dp is the Fourier transform of the test function h{p). 



A comparison of the STF (1.71) with the trace formula (|1.4g|) reveals that both relations are almost 



identical. If one were to set up the trace formula ( 1.48| ) for the systems presently under study one 
would first of all notice that due to the STF all higher-order corrections to the periodic orbit sum 
which are not explicitly contained in ( |1.48D indeed vanish. One is then furthermore able to read off 
from the STF a number of quantities entering the semiclassical trace formula: 

1. The mean spectral density is completely known, 

A A 

do{p) = ttP tanh(7rp) ~ 7— p + 0(6"^"^) , p ^ 00 . 

zvr zvr 



Since the surface has no boundary, in a comparison with ( 1.43 ) one has to choose L = 0. 

2. The eigenvalues of the monodromy matrices M-^ are given by e^"^ = e^'^, that is u^ = l^, for 
all periodic orbits 7. 

3. No Maslov phases occur, e~*2^7 = 1. 

4. The width of the strip of holomorphy demanded for h{p) is explicitly known, and hence P{2) = 2- 
Indeed, a further study of the classical dynamics yields that the topological pressure is linear, 
P{P) = 1-/3. 
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Of course, all of the applications discussed for the semiclassical trace formula in section ^^ carry over 
to the STF. 

All the information on the distribution of classical periodic orbits and on the quantum energies 
which is contained in the STF can be encoded in a function of the complex variable s = 2 — ip, where 

p = \ E — T is the momentum variable introduced above, but now extended to the whole complex 



plane, see |23, 20|. This Selberg zeta function 



Z{s) := n n (^ " e-("+")'^p 

7p 71=0 



(1.72) 



is defined by a product over the classical primitve periodic orbits 7p which converges in the right 
half-plane Res > 1. Since the latter condition is equivalent to Imp > ^ the domain of convergence 
excludes the real momentum axis where all quantum energies, with the exception of possibly finitely 
many, are to be found as En = p1^ + \- 

In order to obtain the analytic properties of the Selberg zeta function one chooses the test function 



h{p) 
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1 



p2 + (s - i)2 



p^ + {a- \f 



(1.73) 



with Res, Recr > 1, and inserts this into the STF. The Fourier transform of h{p) reads g{u) 



2s- 1 
evaluated, 
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2'' so that the sum over periodic orbits on the r.h.s. of (1-71) can be 
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(1.74) 



2s - 1 ds 
Extracting the contribution of the lowest eigenvalue £"0 = on the l.h.s. of the STF then yields 
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n=l 
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En + S{S - 1) En + (J{(T - 1) 
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1 Z'{s) 1 Z'{a 
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(1.75) 



where iIj{z) = ^logr(z). Up to now we have restricted the above expressions to the domain Res, 
Reo" > 1 in order to keep the integrals and sums convergent. The relation ( 1.75 ), however, allows to 
analytically continue Z{s) to all s G C. To this end we perform the limit a ^ 1 and arrive at 



E 

n=l 



1 1 

En + S(S - 1) ~ ^ 



-7A 



s{s - 1) 27r 



A ^, ^ 1 Z'{s) 



2s - 1 Zis) 



(1.76) 



where 7a denotes an appropriate constant, see [24|. An integration of this expression with respect to 
s results in the product representation |24| 



Z(s) = Z'(l)s(s-l)e 



s(s-l)7A 



(2vr) 



[l-s) ^-s(s-l) 



G{s)G{s + l] 



2ir 



n 

n=l 



sis - 1) \ «("-!) 



En 



(1.77) 
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for the Selberg zeta function, where Barnes' double F-function 



G{z + l) = {2TT)^e"^'^'' 



n 

n=l 



Z\^ 

1 + - e 



2n 



(1.78) 



appears. The analytic continuation of the Selberg zeta function based on ( |1.75 ) reveals that Z{^s) is 
an entire holomorphic function with zeros at 

• s = ^ ± i-pn of multiplicity dn-, when En = p^ + |, is an eigenvalue of —A of multiplicity d^, 

• s = 1 of multiplicity one, 

• s = of multiplicity 2;^ + 1, 

• s = -k, /c G N, of multiplicity ^(/c + 1). 



We remark that by topological reasons the area of a compact surface of constant negative curvature 
is such that ^ is an integer. The analytic structure of Z{s) can also be directly read off from the 
product representation ( 1.77] ). 



If one substitutes s 1— > 1 — s in ( 1.77 ) and subtracts the resulting equation from (|1.771 ) one obtains 
the functional equation 



Z(l - s) = Z{s) exp <^ -A 



-I 



u tan(7rn) du 



(1.79) 



upon integrating the difference with respect to s. We now introduce a mean spectral staircase 



Nip) :-- 



P ~ A fP 

do{q) dq = -^^ I q tanh(7rg) dq , 



'0 27r 70 

and then evaluate the functional equation ( 1.7S| ) ioi s = ^ — ip, p G 



Z[- + ip 



^ \ o ~ ^p] ^^P { ~27ri N{p)^ . 



Thus the function 



C{p) := Z i--ipj ex.p{-iTTN{p)} 



-'5 



ip 



exp 



-ITT 



N{p) argZ(- -ip) 



(1.80) 



(1.81) 



(1.82) 



is real when p € M and satisfies ^{—p) = S,{p). The zeros of Z{s) at s^ = ^ ±Pn appear as zeros of ^(p) 
at ±Pn- Thus, upon increasing p starting at p = 0, one successively passes through all p„'s related to 
the quantum energies En = pii + \- Hence, the spectral staircase N{p) counts all zeros of the function 
^{p') in the interval Q < p' < p. When passing through pn the sign of £,{p) changes according to the 
multiplicity dn of the eigenvalue En, i.e., ^(p) is multiplied by (— 1) ". A comparison with the r.h.s. 
of ( 1.821) therefore shows that 



N{p) = N{p) H avgZ i- +ip 



(1.83) 



see also |2^ for a detailed discussion of this relation. The r.h.s. of ( |1.83| ) can be interpreted as a 
decomposition of the spectral staircase into a mean part and a part Nfi{p) = — argZ(2 +ip) describing 
the fluctuations of the staircase. 

At this stage we remark that since the zeros pn of the function S,{p) are related to the quantum 
energies, one can introduce the quantization condition 



' " .2 1 



e(p) = ^ p' + -e{Eo,Ei,E2,...} . 



(1.84) 
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In 1 25, 2(:] methods to evaluate for real p the zeta function, or rather its semiclassical analogue for 
billards, which use cassical periodic orbits are devised, and then zeros of Re^(p) are calculated. This 
procedure leads to good approximations to quantum energies and constitutes an alternative to the 



quantization rules presented in section 1.2 



In view of the previous discussion one notices that the Selberg zeta function Z{s) enjoys many 
properties which are similar to those of the Riemann zeta function 

oo ^ 

^('y-=Y.^ = Ili^-P'T' ^ Res>l, (1.85) 

n=l p 

where the product on the r.h.s. extends over all primes p, see [27| for a detailed discussion, ({s) can 
be extended to a meromorphic function with one simple pole at s = 1 and zeros on the real axis at 
s = —2, —4, —6, .... The infinitely many complex zeros are located in the strip < Re s < 1; according 
to the famous Riemann hypothesis these are of the form s„ = 2 i^^m ^n > 0. One can hence introduce 
the counting function 

Nr{T) := #{n; < t„ < T} = Nr{T) + ^argC (^ + ^^) > (1-86) 

where Nr(T) = 2;^ log 2^~27r~'~^~'~^(T) ^°^ T ^ 00. The Riemann zeta function also satisfies a 
functional equation, 

SO that 

eij(t):=cQ-it) e-2-^«W (1.88) 

is real and even for t E E, ^_r(— t) = Cij(^)- 

The close similarity between the Riemann zeta function (^(s) and the Selberg zeta function Z(s), 
compare in particular ( 1.83| ) and ( 1.86| ), suggests to view the complex zeros s™ = | ±iin as analogues 



of the spectral zeros s„ = ^ ± ipn of Z{s). Thus the quantities t^ + | might serve as 'model quantum 
energies' and their distribution might imitate the distribution of the quantum energies En = Pn-^ j- 

2 Eigenvalue Statistics and Periodic Orbit Theory 

As already mentioned in the Introduction, one of the basic assumptions in the field of quantum chaos 
is that one should be able to observe fingerprints of classical chaos in the distribution of quantum 
energies. That is, given the sequence 

0<Ei<E2<E3<... (2.1) 

of eigenvalues of some quantum Hamiltonian H, the distribution of the quantum energies is expected 
to follow certain universal rules which allow to identify whether the corresponding classical dynamics 
are chaotic or regular. In order to investigate statistical properties of a spectrum suitable quantities 
that reflect the distribution of eigenvalues are required. Obviously, the spectral staircase 

N{E) := # {n; En < E} = / d{E') dE' (2.2) 

JO 

encodes all information on spectral statistics and thus turns out to be a fundamental object for our 
further purposes. Due to the trace formula ( |1.38D for the spectral density d{E), 

N{E)=N{E)+Nfi{E) (2.3) 
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arises as a natural decomposition of the staircase function into a part N(E) that derives from the 
contribution (ll.2l| ) of t = to the trace formula, and a further contribution Nfi{E) reflecting the effect 
of the non-trivial classical periodic orbits. When discussing the Selberg trace formula we noticed that 
for certain classically chaotic systems the contribution of the classical periodic orbits to the spectral 
staircase is contained in the phase of the Selberg zeta function via Nfi{E) = -aigZi^^ + ip) with 
E = p"^ + ^, see ( |1.83| ). Since N{E) essentially gives the volume of that part of the classical phase 
space where H{p,x) < E, see (|1.2| ), the interesting details of the spectral statistics are contained 
in Nfi(E); the latter describes the fluctuations of the staircase about its mean behaviour N{E). It 
therefore proves useful to defold the spectrum, i.e., to introduce the scaled eigenvalues Xn '■= N{En)- 
By abuse of notation we now denote the counting function for the x^'s by 

N{x):=H^{n; Xn<x} = x + Nfi{x) . (2.4) 

The corresponding spectral density then reads 

d{x) =^5{x-Xn) = l + dfi {x) , (2.5) 

n 

SO that the scaled eigenvalues x„ have a unit mean separation. Defolded spectra can thus be directly 
compared in their statistical properties. 

2.1 Measures for the Distribution of Eigenvalues 

The basic idea behind spectral statistics is to view the list xi,X2,x^, . . . of scaled eigenvalues of a 
quantum Hamiltonian as a sequence of random events. The spectral staircase function ( p.2p is then 
considered as a sample function of a random process. The ultimate goal therefore is to characterize 
and classify all possible such random processes in terms of properties of the -classical and/or quantum 
mechanical- physical systems giving rise to the defolded spectra. 

At first one might wonder what can be random about a spectrum of a fixed quantum Hamiltonian? 
In order to answer this question one should recall the mathematical concept of probability. In general, 
the starting point is a probability space, i.e., a measure space with a normalized measure. Any 
measurable real valued function on a probability space is then considered as a random variable. 
Now, N{x) obviously is an integer valued and piecewise continuous function on the positive real line 
R_l_ = (0, cxd). In order to turn the latter into a probability space one needs to define a normalized 
measure. Since M+ is unbounded one cannot simply use Lebesgue measure because it is not finite. 
Instead, one restricts attention to a bounded region on the positive real line, say to an interval 
Ix '■= [x — Ax, X + Ax] of width 2 Ax that is centered at x. In order to obtain information on the 
asymptotic distribution of eigenvalues one then has to consider the limit x ^ c«. As this is equivalent 
to the semiclassical limit h ^ the asymptotic analysis of spectral statistics allows for an application 
of semiclassical methods as, e.g., the use of trace formulae. 

In defolded spectra the mean separation of eigenvalues is one so that in the mean the intervals Ix 
contain 2Ax quantum energies. The latter number should approach infinity in order to yield reasonable 
limit distributions as x — > oo. We therefore choose Ax = ex", with some suitable constants c > and 
a. If not stated otherwise, a will be fixed to obey ^ < a < 1. In view of the limit x — > oo we intend 
to perform, the unboundedness of the spectral staircase suggests an alternative choice for the random 
process to be studied as a means to characterize spectral statistics. One customarily introduces the 
number 

nL{y):=N{y + L)-N{y) , y e Ix , (2.6) 

of eigenvalues in intervals (y, y + L] of length L > 0. On the interval Ix, from which y has to be 
chosen randomly, we make the simplest choice for a probability measure, dfixiv) '■= 2kx^ix{y) d-V^ 
where Xixiv) denotes the characteristic function of the interval Ix- Given an integrable function f{y) 
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on Ix its mean value with respect to dfixiu) is thus given by 



EM] ■■= I f{y) d^y) = ^j2^^ fiy) dy ■ (2-7) 

One can choose as a probabihty measure on I^ any normahzed measure of the form w{y) dy, where 
w{y) > is continuous and vanishes outside Ix- The hmit distributions then induced by the random 
variable nL{y) on J^,., for x -^ oo, do not depend on the choice of w(y), see |28| for a discussion on this 
point. We therefore stick to our favourite choice dfix{y)- 

At this point we introduce some notation that will be kept throughout the rest of this presentation: 
Above we indicate the interval Ix by a subscript x which we attach to the random variables and their 
distributions. When passing to the limit j; ^ oo, the resulting distributions take care of the asymptotic 
properties of a given spectrum. All quantities corresponding to this limit will then carry no subscript. 
The same notation will be used in case a random process is stationary, i.e., no dependence on x occurs 
anyway. 

Now, for each value of the parameter L the function nL{y) is a random variable on the interval Ix, 
and thus upon varying L it defines a continuous parameter random process. The ultimate aim then is 
to investigate the resulting limit random process as x — > cxd as far as possible. On the way, however, 
we will notice that already at finite x one can obtain interesting information. As long as L is finite, 
nL{y) provides information on correlations among approximately L eigenvalues and therefore yields 
the local spectral statistics. In order to reach the global scale one has to perform the limit L -^ oo. 
Starting with nL{y) on Ix one is therefore confronted with a competition of the two limits x — > oo and 
L — >■ oo. Phenomena occurring in this context are the central topic of the third part of these lectures. 

For the further analysis it proves useful to introduce the probability Ex{k; L) to observe a value of 
k for the random variable nii^y) when y is uniformly distributed on the interval Ix, 

fX+Ax 



-I pX + £\X 

Exik;L) := lEx [(5„^,fc] = — r— / ^nUy),k dy , (2. 

ZL\X Jx-Ax 



where dk^i denotes Kronecker's symbol. According to ( |2.5| ) and ( p.7[ ) the mean value of riLiy) reads 
lEx[nL] = ^ /^^ " [^(y + L)- N{y)] dy = L + ^ J^^^^^ [Nfi{y + L) - Nfi{y)] dy . (2.9) 

A simple result emerges in the limit x ^> oo because Nfi(y) describes the fluctuations of the spectral 
staircase about its mean behaviour N{y) = y. Thus in particular 

hm -^ / Nji{y)dy = {), (2.10) 

x^oo 2Ax Jx-Ax 

so that the mean value IEx[nL\ asymptotically approaches L. See [p^ j for a discussion of ( p.lOj ) which 
makes use of the trace formula. Having in mind the limit x — > oo as our actual goal, we introduce the 
m-th moments of the shifted random variable ni{y) — L, 

S-(L) := Ex [{UL - in = ^ £^^' my + L) - N{y) - L]^ dy . (2.11) 

Due to the above remark S;|,(L) vanishes as x — > oo so that the lowest non-trivial moment is given by 
the so-called number variance 

^'^^^ = 2A^ / ^^^y + ^^ " ^^y^ ~ ^1' ^^ = ^^^ " ^^' ^"^^' ^^ • ^^-^^^ 



fc=0 
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One can easily determine the behaviour of the number variance for small L: Since Ix contains finitely 
many scaled eigenvalues there exists a positive minimal separation s„ := x„+i — x„ of neighbouring 
levels. Once L is smaller than the minimal Sn every interval {y, y + L] with y € Ix either contains one 
or no Xn, and thus Ex{k; L) = for k > 2. One then employs the normalization of the probabilities, 
ExiO; L) + Ex{l; L) = 1, to observe 

1 
S2(L) = Y,{k - L? Ex{k; L) = (1 - 2L) iE,K] + L^ . (2.13) 

fc=0 

The relation ( |2.9| ) then implies that 

T,l{L) r^ L , L ^ . (2.14) 

The behaviour of 5]^.(L) for small L does not provide detailed information on the distribution of 
eigenvalues since for too small L the number variance merely reflects the fact that N(x) is a staircase 
function. 

The opposite limit L — > 00 cannot so easily be dealt with, but on the other hand reflects specific 
properties of the underlying physical system. From ( ^.12 ) we obtain 

^1{L) = -^ / [N{y + L)-N{y)-Lfdy 

= Ex+L [Nji] + Ex [N]i\ - 2Ex [Nfi{y + L)Nfi{y)] . (2.15) 

Certainly, the limit L — > 00 also requires to perform x — > cxd, since otherwise the intervals (y,y + L] 
for y € Ix hardly overlap with Ix] by the same reason L should not grow faster than x. We therefore 
choose L = L[x) = ax"' with some positive constants a and 7 < 1. Now two effects on the behaviour of 
S^(L) in this limit are to be observed on the r.h.s. of ( p.l5| ). On the one hand, a large exponent 7 will 
decouple Nfi{y + L) and Nfi{y) to a certain extent so that the contribution of the correlation function 

on the r.h.s. of ( p.l5[ ) becomes smaller in absolute value. On the other hand, a small exponent 7 will 

-2 
fi 



turn the two expectation values of Nf, almost equal. However, in the asymptotic regime for x — > 00, 



7 need not really be small because once Ex[N'ii] increases for x ^ 00, then lEx+LWfi] ~ ^xWfi] for 
any choice L = ax"' with 7 < 1. It indeed turns out that in typical cases the expectation values of 
N%{y) on Ix behave either like 

Ex [N]i] ~ 61 logx or Ex [N]i] ~ 62 x" , (2.16) 

with some /) > 0. We hence expect that as long as 7 is larger than some critical exponent 7c the number 
variance is asymptotically given by the dominant and L-independent contribution 2Ex[N'ij\, with L- 
dependent modifications provided by the correlation function Ex[Nfi{y + L)Nfi{y)]. Theoretical [^] 



as well as numerical fsy, ^ studies indeed confirm this qualitative picture. In section |2^ we will 
return to a more detailed discussion of this point. 

The characteristic function of a probability distribution is defined as the Fourier transform of the 
probability measure. In the case of the random variable n/,(y) it is given by 



1 fX+Ax 



Jn.AO ■■= JEx e^^"^ = ^A^ y e^'""^'^ ^2/ • (2-17) 



The probabilities Ex{k;L), see (|2.8| ), can be recovered from the characteristic function because 

1 /-x+Ax 1 px+Ax 1 r+TT 

= ^ r'e-'(''Jn,Ai)di- (2.18) 



— 7T 
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Hence the knowledge of the characteristic function is equivalent to knowing the probabilities Ex{k; L) 
themselves. Characteristic functions can be used to test the statistical indepedence of random variables 
since the characteristic function of the sum of independent random variables factorizes into the product 
of the respective characteristic functions of the individual random variables. In our case, ni{y) = 
L + Nfi(y + L) — Nfi{y) essentially is the sum of the two random variables Nfi{y + L) and —Nfi{y), 
with y € Ix- If one wants to test their statistical independence one therefore has to study whether or 
not the characteristic function ( |2.17| ) factorizes. Certainly, a complete statistical indepence can only 
be expected for the limit distributions as x — > oo, since only these take infinitely many eigenvalues 
into account. As argued above, a decoupling of the distributions of Nfi{y + L) and Nfi{y) will then 
require to perform also L ^ oo, with L = ax"' and a positive exponent 7 which is large enough. One 
way to determine the minimal such 7, which then yields 7c, could be to test the factorization of the 
characteristic function Jni^,x{C)- At the end of section |2.3| we will return to this topic in some more 
detail. 

2.2 Random Matrix Theory 

So far we discussed the distribution of eigenvalues of individual quantum Hamiltonians, without speci- 
fying any result about, e.g., the probabilities Ex{k; L). However, we already addressed the expectation 
to observe fingerprints of the classical dynamics on the corresponding distribution of quantum ener- 
gies. Ideally, this could mean a universal behaviour of the probabilities Ex{k; L), as x — > 00 of course, 
depending only on the type of the classical dynamics. For example, one could expect that all classi- 
cally chaotic systems lead to one and the same limit distribution for their respective random variables 
riL', classically integrable systems then should produce a distinctively different distribution. This ex- 
pectation being very strong could be weakened in replacing 'all systems' by 'almost all systems' in 
order to account for 'special' or 'non-generic' cases where exceptions might occur. Of course, one 
then has to specify what 'almost all' means in this context. To this end it is necessary to consider a 
family of quantum Hamiltonians Hx, where A ranges over some parameter space A. Since now one is 
dealing with families {£'i(A), E2{X), E3{X), . . .} of quantum spectra all spectral quantities depend on 
the parameter A, e.g., nL^\{y), Ex^\{k\ L). Provided the parameter space A is endowed with a suitable 
measure di^{X), a weak form of universality in spectral statistics would be reflected in the fact that 
the limit probabilities Ex{k;L), arising from Ex^x{k;L) as x — > 00, are given by Ep^{k;L) for almost 
all A G A with respect to the measure dv[\). If the latter were normalized, j^dv{\) = 1, i.e., it were 
a probability measure, then one would obtain 

Etv{k]L) = I Ep,{k;L)dv{\) = I Ex{k-L)dv{\) . (2.19) 

JK Ja 

Thus, even in case it is not known whether Ex{k;L) = E/<^{k;L) for almost all A G A it would be 
interesting to calculate the mean value on the r.h.s. of ( p.l9|) which describes the spectral statistics of 
the quantum Hamiltonians Hx on average. In this perspective the limit of Ex^xik; L) as x ^ 00 may 
even not exist for a set of A which is of zero measure in A. 

Examples for families of classically chaotic systems to which the above programme might be applied 
are provided by geodesic motions on surfaces of constant negative curvature as discussed in section 
[1.3| . The compact such surfaces with the topology of a sphere with g handles {g > 2) form a {6g — 6)- 
parameter family. For g = 2 Aurich and Steiner ||32[| calculated numerically eigenvalues of quantum 
Hamiltonians on 30 surfaces which they picked out randomly from the corresponding 6-parameter 
family, and computed averages of the probabilities Ex^x{k]L) over their sample of 30 values for the 
parameter A. Apart from the fact that x necessarily had to be kept finite and thus the saturation 



effects to be discussed in section ^^ showed up, it was observed that the 'ensemble averages' agreed 
with the respective probabilities in the Gaussian orthogonal ensemble (GOE) of random matrix theory. 
We remark that in this example it turns out to be essential not to expect the GOE result for all values 



of A since certain exceptions, which go under the notion of arithmetic quantum chaos [33, 34, ^, are 
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known to exist. However, with respect to any reasonable measure dv{X) these will form a set of zero 
measure. 

In a purely mathematical model the programme has recently been completely carried out by 
Katz and Sarnak [ p7[ . As model operators they consider matrices from the classical compact Lie 
groups 0{N), SO{N), U{N), SU{N), and USP{N). Then the respective group itself can serve as 
the parameter space A, and the most natural choice for dv{\) is provided by the Haar measure on 
the group. Katz and Sarnak prove that in the limit A^ — > oo almost all matrices from any one of 
the classical groups show the spectral statistics of the Gaussian unitary ensemble (GUE) of random 
matrix theory. 



In random matrix theory, see [38| and Bohigas' contribution to these proceedings, one attempts 
to construct from certain assumptions the probabilities E\{k;L) which are expected for 'almost all' 
quantum Hamiltonians of a given type. In this context the first step consists of cutting off the 
dimension of the Hilbert space on which a Hamiltonian H is defined at a finite value A^. Then, upon 
choosing some orthonormal basis (onb) in the Hilbert space, H is represented by a hermitian A^ x A^- 
matrix M. In case the physical system from which IH derives is time-reversal invariant one can always 
choose the onb in such a way that IH is real symmetric. For simplicity we will in the following always 
focus our attention to this case. One then still has the freedom to choose another onb without changing 
H and its spectral properties; only the matrix representation is affected. The base changes in the A^- 
dimensional Hilbert space are provided by orthogonal transformations O € 0{N). The matrix M is 
thus conjugated to 0~^IHO. As a model space of truncated Hamiltonians one therefore considers the 
set of all real symmetric matrices M, with an identification of M and O^^MO whenever O € 0{N). 
Altogether a real symmetric A^ x A^-matrix has iA^(A^ + 1) independent entries. The symmetry group 
0{N) is of dimension ^N{N — 1) so that the set of classes of equivalent real symmetric matrices has 
dimension A^. A possible choice for coordinates of the space of truncated Hamiltonians is obtained 
after diagonalizing each M: There exists an O G 0{N) such that 0~^1H0 is diagonal. Therefore 
each IH is equivalent to a diagonal matrix, with its eigenvalues Ei, . . . , Ejsf as entries on the diagonal. 
The eigenvalues then yield A^ parameters to describe M, and the additional ^^{^ ~ 1) coordinates 
required to characterize M completely are provided by the independent entries of O. 

The basic assumption of random matrix theory now concerns the measure dvi^{]H) on the space 
of real symmetric A^ x A^-matrices. It is based on the principle of minimal knowledge and consists 
of two parts: (i) dvN{JH) = PN{IH)dlH, where dH = UkKi^^^ki, and Pn{IH) = Pn{0-^HO), so 
that the measure does not depend on a particular choice of an onb. (ii) The independent entries M^i^ 
k < I, are statistically independent. By these requirements the measure dviq{]H) is already fixed to 
read |8| 

dvN{JH) = const. eyi.-p {-ati M'^ + hii H] dH , (2.20) 

with some real constants a > and h. In order the obtain a more explicit expression one changes 
variables from M^u k < I, to Ei, . . . ,E]\f and 2N{N — 1) further variables which parametrize the 
orthogonal matrix O that diagonalizes M. One finally observes |38] 



N 

?2 



duNiiH) = dp{o) n i^fc - ^'1 n ^~"^' ^^i ' (2.21) 

k<l j=l 

where dp{0) denotes some measure that only depends on O. Upon integrating out dp{0) one has 
thus obtained a measure on the model space of truncated Hamiltonians, where the latter are being 
parametrized by their eigenvalues. Notice that the eigenvalues are not statistically independent since 
the measure does not factorize. The term nfc<zl-^fe ~ Ei\ rather introduces a correlation of the 
eigenvalues that makes small separations E^+i — E^ of neighbouring eigenvalues less probable. This 
effect is known as level repulsion. 

The limit measure resulting from ( p. 21 ) as A^ — > cxd is commonly called Gaussian orthogonal 



ensemble (GOE). Given the latter one can then calculate all statistical properties of eigenvalues of 
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random symmetric matrices. For example, the number variance reads 



^goe{L) = ^ [log(2^L) + 7 + 1 - cos(2^L) - Ci(27rL)] 



vr- 
+2L 

2 



vr 



-Si(27rL) 



1 



1 



+ ^Si(7rL)^ - -Si(7rL) 



vr^ 



vr 



vr 



(2.22) 



log(2vrL) + 7 + 1 - ^ 



.0,1 



where 



Si(x) := 



smy 



dy and Ci(x) := 7 + logx + 
y JO 



cos y — 1 



dy 



(2.23) 



Here 7 = 0.577215 . . . denotes Euler's constant. For the GOE one does not know, however, an exphcit 
closed expression for the probabilities EcoEik; L) to find k defolded eigenvalues in spectral intervals 
of length L; but these quantities can be calculated numerically to a high precision. That way Aurich 



and Steiner [g2| found 



(k-Ly 



EGOE{k;L) 



27:TIoe{L) 



"GOE 



(L) 



L 



00 



(2.24) 



Later Costin and Lebowitz JS^] made precise the way the limit of large L has to be understood. 
Since both the mean value JEcoEinL] = L and the variance ( |2.22| ) of the random variable n/,, when 
considered in the GOE, diverge as L — > 00 one has to employ a renormalization before one performs the 
limit. Costin and Lebowitz proved that after renormalization the distribution of the random variable 



obtained from n^ weakly approaches a standard Gaussian, see section 3.1 for a further discussion. 

The construction of the GOE relies on a number of assumptions that cannot be derived from 
first principles when one is discussing the distribution of eigenvalues of actual quantum Hamiltonians. 
Historically, random matrix theory was introduced to describe the statistics of energy levels and of 
resonances of large atomic nuclei. In view of the huge number of degrees of freedom involved the 
principle of minimal knowledge seemed to be reasonable in this context. It indeed turned out to result 
in a correct description of experimental data, see |40| for a collection of the original contributions. Later 
Bohigas, Giannoni, and Schmit ||4^ conjectured that the eigenvalue statistics of quantum Hamiltonians 
with even a low number (> 2) of degrees of freedom universally follow the predictions of random 
matrix theory, if only the corresponding classical dynamics were chaotic. Again, the conjecture found 
confirmation by numerical calculations of eigenvalues and their statistics in many examples. 

For classically integrable systems the situation appears to be completely different. Since at least 
a second classical integral of motion exists one has a quantum mechanical observable at hand that 
commutes with the Hamiltonian. For each value of the quantum number that corresponds to the 
second conserved quantity one obtains a subspectrum of the quantum energy spectrum. It now seems 
that the entire energy spectrum behaves like an independent superposition of these, typically infinitely 
many, subspectra. Provided the superposition of subspectra were indeed uncorrelated, the complete 
spectrum would have the local statistics of a Poissonian random process, i.e.. 
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(2.25) 



In particular, the density P{s) of the distribution of nearest neighbour level spacings Sn := Xn+i — Xn 
of defolded eigenvalues would be exponential, Ppoisson{s) = e~*. Based on a semiclassical analysis 
the Poissonian behaviour of the local eigenvalue statistics was predicted by Berry and Tabor [^ and 
meanwhile has been verified by numerical calculations in many examples. See e.g. P3| for a careful 
numerical analysis. 
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But again, exceptions are known to exist. In order to give an example suppose that ei,e2 S M? 
are linearly independent such that C = Zei + Ze2 is a lattice in M?. Then ]R^/£ is a two dimensional 
flat torus. The classical geodesic motion on M.'^/C is integrable, and the spectrum of H = —-^A on 
L^(M^/£) is explicitly known; it comprises of the eigenvalues Ep;^ = 2^4:Tr'^{kfi + 1/2)"^- Here k and I 
run through Z, and {/i, /2} is a basis for the dual lattice £* of C, i.e., Cj • /,• = (5jj. Hence 

where A,B,C are suitable real constants which determine Q{k,l) to be a positive definite binary 
quadratic form. If -^ and ^ are rational numbers, Q{k, I) is a rational quadratic form. In this case the 
value of Q{k, I) itself is a rational number for all k,l £ Z; it can in particular be an integer, Q{k, I) = n, 
n G N. If we denote the number of pairs {k,l) G Z^ that represent n as n = Q{k,l) by rQ{n), the 
multiplicity of the eigenvalue E = 27^ '^ is rQ{n). Now, for rational forms it is known that rQ{n) is 
unbounded as n — > 00. One can then conclude that the distribution of level spacings becomes singular 
at s = 0; it is actually given by the 'density' P{s) = 6{s), see for example [p^]. 



The coefficients ^, i?, C in ( 2.26 ) provide parameters for the set A of all quantum Hamiltonians H\ 
that derive from Laplacians on flat tori. The positivity of the eigenvalues only requires the discriminant 
Dq = B^ — AAC to be negative. The set of Hamiltonians in A corresponding to rational quadratic 
forms then clearly is of zero Lebesgue measure, so that the tori leading to the singular distribution 
of level spacings as described above are non-generic in a well defined sense. It is widely believed, 
however not known, that for a set of full Lebesgue measure in A the level spacings distribution has an 
exponential density, P{s) = e~*. As an important step towards a proof of Poissonian local eigenvalue 
statistics for Laplacians on flat tori Sarnak recently proved |^^ that the spectral statistics, when 
measured with the pair correlation function, is Poissonian on a set of full Lebesgue measure in A. 
Moreover, if i? = and ^ satisfies a certain Diophantine condition, Bleher and Lebowitz [45] proved 
that the number variance S^(L) is Poissonian. In view of the above discussion a further result of 
Sarnak is rather amazing. He could show that there exists a set A C A of tori, which is generic in 
a topological sense, and on which the level spacings distribution has no exponential density; in fact 

-1/0 ~ ~ 

Jq Px{s) ds does even not exist for A S A. The set A is generic in the sense that it is of second Baire 
category, i.e., it is a countable intersection of open dense sets. The lesson to learn from this example is 
that universality in spectral statistics can only hold in a weak form, as for example expressed around 
( |2.19| ). One even has to be very careful in the way one is declaring particular quantum Hamiltonians 
as 'generic'. 

2.3 A Semiclassical Analysis of Spectral Statistics 

It seems that the only tools available to investigate the distribution of eigenvalues of a single quantum 
Hamiltonian H are provided by semiclassical trace formulae. These relate the spectrum of H to classi- 
cal quantities: periodic orbits and their stabilities and Maslov phases. One therefore might anticipate 
that the asymptotic distribution of eigenvalues is completely determined by the corresponding classical 
dynamics. The first to perform such a semiclassical analysis, to classically integrable systems, were 



Berry and Tabor [42|. Later, Berry extended this analysis to classically chaotic systems as well, see 
[^ , Q and his contribution to j^. Here we only want to consider quantum systems with a strongly 
chaotic classical limit, which shall in particular mean that all classical periodic orbits be unstable 
and isolated in phase space. To be definite, we further restrict our attention to cases where the trace 
formula ( 1.4§| ) applies, and moreover, the spectral density has an asymptotic expansion as in ( |l.43| ). 



Thus, the classical dynamics are given by either billiards or geodesic motions on surfaces. The case 
of the Selberg trace formula (1.71|) is therefore included. 



25 



The starting point for our discussion is the representation 

r-y+L ry+L 

nL{y) = N{y + L)-N{y)= d{u) du = L + dfi{u) du (2.27) 

Jy Jy 



of the random variable nL{y), y £ Ix, in terms of dfi{u), compare ( p.SI ). Now, the fluctuating part of 
the spectral density, when expressed in terms of the momentum variable p, can be given in terms of 
a periodic orbit sum as in ( [1.59 ). The number variance, being the second moment of ^^(y) — L, can 
hence be represented as 



fX+Aa; 



S2(L) = ^r ^[nL{y)-Lfdy 

1 /-x+Ax ry+L ry+L 

— — / / / dfi{ui)dfi{u2) duidu2 dy. (2.28) 

^^ Jx~Ax Jv J y 



Zl\X Jx~Ax Jy Jy 

To proceed further, one customarily introduces the two-level form factor 

hoo 1 /-x+Ax 



g2.itr / My) My + 1) ^y at. (2.29) 

-00 ZZAX Jx-Ax 

By Fourier inversion and a change of variables one obtains from ( p.28| ) and ( 2.29 ) that 

/+CXD pL pL 

dr dpi dp^e-^-^^P--"-^ Kx+p,{t) . (2.30) 

-oo Jo Jo 

As we are finally interested in the limit of x ^ oo, with < pi < L <^ x, we approximate /fa;+p^(r) ~ 
Kx{t) and observe 

E^(L).4r!i!^A-(r)d.. (2.31) 



VT" Jo T- 

After having performed re — > oo the above asymptotic relation indeed turns into an identity. When 
differentiating the latter with respect to L and applying a sine-Fourier inversion afterwards, one obtains 



K(r) = ±^HL) 



L=0 



1 f°° d^ 

— / sm(2TT Lt)— -,J:^(L) dL , (2.32) 

TTT Jq dL'^ 



from which one concludes that 



K{t) = 1 + o(-j , r^oo, (2.33) 

since ^^(L) ~ L for L ^ 0, see (|1|). 

We are now going to express the form factor in terms of sums over classical periodic orbits. To this 
end dfi{x) in ( |2.29| ) has to be represented by the periodic orbit sum (|l.59D . However, since the form 
factor ( |2.29| ) is given in terms of the defolded spectral variable x = N{E), we first have to rewrite 
it as a function of the energy E = p^. From the relation ( |l.43| ) we conclude that N{E) ~ j:;^^ ^^ 
S — > oo or ?i — > 0, and thus d{E) ~ j-jrj =: d. We therefore observe that 

, / X d ^^ , , dE dNfi(E) dfi(E) --i , ,^, , ^, 

The two- level form factor, when expressed by E instead of x, hence reads 



p+oo _ 1 i-E+AE / >.\ / y 

K{t- E) ~ d~' I e^--'^^ -— / dfAE'--]dfAE' + -] dE' dX , (2.35) 
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where we now correlate dfi at the symmetrically placed points E' ± -^ instead of E' and E' + A. 

We are now in a position to employ the semiclassical periodic orbit sum ( |1.59| ) for dfi{p) 
2pdfi{E), p = vE, in ( 2.35 ). The result, which was first obtained by Berry |^^, then reads 



K{t;E) 
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where 
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(2.36) 



(2.37) 



2P |det(M;^-l)|2 2p 

is an amplitude attached to the r-fold repetition of the primitive length In, which is of multiplicity 
gn- The inner double sum in ( |2.36D extends over all non-zero integers r, s, both positive and negative. 
The braces < . . . > appearing in ( p. 361) denote the average in E' over the interval [E — AE, E + AE]. 
Furthermore, the quantity Ih '■= 4:Trphd is introduced. Apart from the factor of 47r, which is subject to 
convention, Ih is a combination of quantities which has the dimension of a length and can be used as a 
'semiclassical parameter' in that the semiclassical limit corresponds to Ih — >■ oo. (The subscript refers 
to the fact that Ih is sometimes called 'Heisenberg length'.) Notice that in the discussion following 
( |1.63| ) Ih has been identified as the length where the periodic orbit sum ( 1.63| ) can be truncated without 
leaving out its most important contributions. Following this philosophy, we therefore also truncate 
the periodic orbit sums in ( |2.36| ) at |r/„|, \slm\ < Ih- Then obviously the resulting expression can at 
most represent the form factor in the range \t\ < 1; for larger values of |t| the truncated periodic orbit 
sum yields no contributions. But already in the interval ^ < t < 1 some terms are missing so that a 
reliable statement can in fact only be made for |r| < ^. Since now we are dealing with finite sums the 
semiclassical representation of K(t; E) is no longer plagued with convergence problems. The price to 
pay for this convenient fact is the restriction to |r| < ^, and the inaccuracies which derive from the 
omission of the -infinite- tails of the periodic orbit sums. However, now the non-diagonal terms on the 
r.h.s. of ( p. 36 ), with r/„ ^ slm, are suppressed by the average < . . . > over E', or equivalently over 



p . In the limit E, AE -^ co the non-diagonal terms actually vanish. We have therefore arrived at the 
semiclassical representation 

K{t;E) 
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(2.38) 



for < r < i and E ^ oo. The subscript refers to the fact that Kdiag{T]E) is solely defined by 
the diagonal terms of the periodic orbit sums. Due to the exponential proliferation of the number of 
periodic orbits 7 with lengths l^ < /, see (1.64), the sum over r > 1 is dominated by its contribution 
from r = 1. If we moreover substitute for the amplitude An^i its leading asymptotic behaviour ( |2.37| ) 
we obtain 

Kir; E)r.^y^glll e-«" slr-^-^^ . (2.39) 
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The sum over the length spectrum in ( 2.39 ) will now be rewritten as a sum over periodic orbits, with 
the effect that one factor of Qn is being absorbed by the sum. We furthermore now assume that the 
multiplicities Qn of lengths /„ asymptotically approach a constant value of 'g. Thus 



K{r;E)^^ Y: 



l',e- 
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(2.40) 



We then integrate K{t';E) with respect to r' G [0, r], r < ^, 



K{t', E) dr' 
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(2.41) 
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Periodic orbit sums like the one appearing in (2.41) frequently occur in the thermodynamic formalism 
for hyperbolic dynamical systems, see for example |^3l- The asymptotic behaviour as /// — > oo of the 



r.h.s. of (2.41) can be obtained from relations known in the thermodynamic formalism, which relate 
the various dynamical entropies and exploit properties of the topological pressure -P(/3), see ( |1.51| ). 
For hyperbolic systems with compact phase space one can in particular show that P{1) = 0, and from 
this result one can conclude that ( |2.41 ) yields |t^, so that after differentiation with respect to r, 



K{T-E)^-gT 



T< 



E 



oo 



(2.42) 



The linear behaviour of the form factor was first obtained by Berry [^]. He concluded the result 
( |2.42| ) from ( 2.40| ) with the help of a sum rule due to Hannay and Ozorio De Almeida |48|. 

Above we derived the general behaviour ( 2.33 ) of the form factor as r -^ oo. Since it must 
approach one, the linear increase ( 2.42| ) cannot extend beyond r ~ 2- We recall that the restriction 
to r < 2 resulted from the truncation of the periodic orbit sums at Ih- It now becomes obvious 
that for large r the r.h.s. of ( p. 36 ) has to behave in a way that cannot be controlled by as simple 
semiclassical considerations as for small r. A further observation that can be made with the relation 
( 2.42| ) is that the behaviour of K{t;E) depends on the average asymptotic multiplicities of lengths 
of periodic orbits. In 'generic' classical dynamical systems with time-reversal invariance one obtains 
that g = 2 because a periodic orbit and its time-reversed image share the same length. Here 'generic' 
means, among other things, that the classical dynamics have been completely desymmetrized, i.e., all 
discrete symmetries have been removed; a multiplicity of (7„ > 3 is then of an accidental kind. And, 
moreover, only if an orbit of length In is self-retracing, one observes gn = 1. Typically, the asymptotic 
fraction of self-retracing orbits vanishes so that indeed 'generically' 'g = 2. For small /„, however, 
effects of multiplicities g„ = 1 remain. Hence, the relation ( 2.42| ) can only apply once one is deep 
enough in the semiclassical regime, i.e., for sufficiently large Ih- When the classical dynamics shows 
no time-reversal invariance, one typically finds g„ = 1 so that then g = I. Thus quantum systems with 
'generic' time-reversal invariant classical limits will yield a number variance that differs from those 
quantum systems whose classical limits are not time-reversal invariant. In random matrix theory a 
lack of time-reversal invariance forces one to deal with complex hermitian matrices and unitary base 
canges. Then the measure replacing the GOE is the Gaussian unitary ensemble (GUE). 

At this point we can also hint at the reason for the exceptional behaviour of the spectral statistics 
in arithmetic quantum chaos, see 135 
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, IS^l- There the classical systems are completely chaotic, 
but the multiplicities of lengths of periodic orbits are exceptionally large; indeed one can show that 



5n ~ c- 



I'll 



for L 
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]. Clearly, the above discussion of the form factor has to be modified 
in an essential way. As a result, one can show [29| that K{t;E) increases exponentially instead of 
linearly. 

We are now going to analyse the properties of the number variance by means of the relation ( |2.31| ). 
The form factor will be modelled by a simplified version that, however, captures its main features 
as they have been worked out above. For small r, close to zero, the form factor will be represented 
by its diagonal approximation (2.38), which we deliberately will sometimes use in the form ( |2.39 ) or 
( 2.401 ). For larger values of r, however still below 2, K{t]E) can be well approximated by its linear 
overall increase ( |2.42| ) with g = 2. Then, at r = ^, both the approximations employed to arrive at 
( |2.42| ) break down and the form factor reaches its asymptotic value of one, see ( 2.33| ). For simplicity, 
and because no further information on the behaviour of the form factor beyond r = 2 is available, we 
merely set K{t; E) = 1 in our model. This therefore reads in explicit terms 



Km{t-E):-- 




< r < T* 

T* <T<1 



(2.43) 



see also [^6[. At this point we remark that recently Bogomolny and Keating |51] managed to take 
non-diagonal contributions to ( 2.36| ) into account and hence provided a further justification for ( |2.43| ). 
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The value of t* , beyond which we replace KdiagiT; E) by its semiclassical asymptotics (2.42), seems to 
be somewhat arbitrary. Its actual role will become clearer in the course of our subsequent discussion. 
As a first illustration of the meaning of r* we remark that ( 2.42| ) results from the estimate 



■y,l-y<l 
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oo 



(2.44) 



applied to ( 2.41 ) with / = Iht. As long as the cut-off / is not connected to Ih no quantum mechanical or 
semiclassical quantities appear; ( p. 44 ) thus is a purely classical asymptotics. One can hence introduce 
a classical length scale Id such that the r.h.s. of ( p. 44 ) approximates the l.h.s. to a given precision. 
We then define r* := Id/ln- Thus, for t > t* one obtains that ///r > Id which allows to use the 
asymptotics (2.44). The arbitrariness in the choice of r* in ( p. 43 ) now appears as the -restricted- 
freedom to fix a value for Id- Certainly, Id has to be considerably larger than the length ^i of the 
shortest periodic orbit. Moreover, since t* < -^ and indeed the semiclassical limit corresponds to 
Ih -^ oo, we have three well separated length scales ^i <C Id "^ Ih- 
The model number variance now follows from ( p. 431 ) through 



^lliL;E):-- 
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and consists of three parts which are due to the three r-ranges as they appear in ( |2.43| ). The first 
part derives from the integration over the interval < r < r* and reads 
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(2.46) 



It depends on r* throught the cut-off Id- The second part, which emerges from the integration over 



r* < r < i, is given by 



^M,2 
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[- log(2r*) - Ci(7rL) + Ci(27rLr*)] , 



recall ( |2.23| ) for the definition of Ci(x). Finally, the remaining part reads 
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(2.48) 



An immediate observation that can be made with ( |2.46| )-( 2.4^ ) is that Tj\.^(L;E) ~ L for L — > 0, 
which is in accordance with the general statement ( |2.14[) . 

Having explicit expressions at hand for the number variance, we are now in a position to analyse 
its behaviour for large L. To this end we first notice that S|^(L;£') can be split into a sum of two 
contributions, one of which is explicitly independent of r* or Id-, respectively, and reads 



^M,a(-^;-^) 



The second part 
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(2.50) 
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contains r* and Id- We will see below, however, that the leading order asymptotics as L ^ oo 
is independent of the parameter Id- Inspecting ( |2.5[l| ) one observes that L and r* appear in the 
combination Lt* . Now, r* = Ici/Ih — > in the semiclassical limit Ih — > oo, which forces us to analyse 
the limit L ^ oo in the two cases Lt* -^ and Lt* — > oo. For convenience we now introduce 
the dimensionless quantity L^ax := 7* so that in the following we discuss (i) L/L^ax —>■ and (ii) 

L/Lfnax ^ 00. 

In the first case the argument of sin^ can be estimated as -t^ttL < ttt <C 1 so that the periodic 

orbit sum in ( 2.50 ) yields a contribution of 0{{L/Lmax)'^)- Altogether, one obtains for ( 2.49|) and 



( |2.5[1| ) in the limit L -^ 00, Ih — > cxd, L/Lmax -^ 0, 



-|[log(27rL) + l+7-log2] + On-^j ] + O (\ 



Si,(L;i5;) = — [log(27rL) + l+7-log2]+0( (7 ) )+o(-). (2.51) 



Our first observation with ( |2.51| ) is that, as announced before, all terms that have been evaluated 
explicitly do not depend on l^i- Furthermore, a comparison with the number variance ( 2.22| ) in the 



GOE reveals that the leading terms coincide. Concerning the next-to-leading order, however, the 



^•2 



constant in (2.51) contains —log 2 ~ —0.693 instead of — ^ ~ —1.234 as in the GOE. Thus, our 
model ( ^.43 ) reproduces the number variance obtained in random matrix theory to leading order. 
The next-to-leading terms, however, differ slightly. This effect is most probably due to the crude 
approximation ( 2.43[) to the actual form factor in the domain r > 2- One might expect that the true 



form factor will fully reproduce the GOE result. 

Remarkable deviations from random matrix theory show up in the second case L/Lmax — *■ 00 ■ 
Here the arguments of the sin^-terms are no longer necessarily small; hence the periodic orbit sum 
contributes oscillatory terms. We are interested in the overall behaviour of the number variance as 
L -^ 00, Ih -^ c)o, L/Lmax ^ 00 so that we now determine the mean value Ti^{E) about which the 
oscillations of the number variance take place. This can be achieved by replacing each sin^-term by 
its average value g- Thus the periodic orbit sum no longer contains L; it merely yields a constant 
C{lci) that is solely determined by the value of the cut-off length Id- Collecting then the asymptotics 
of ( 2.491) and ( 2.50| ) in the limit L ^ 00, //f — > 00, L/Lmax -^ co yields 

^ [log Lmax + 1 - log 2 + C{lcl)] + O ( —^ 



S^(^) = — [\0g Lmax + 1 - log 2 + Cilci)] + O ( -^ ] . (2.52) 



It now seems that the result ( p.52| ) depends on l^i- However, the following estimate shows that to 
leading order this dependence disappears. To this end take only the leading contribution of the 
primitive periodic orbits into account, 

Cild) « A E 5"^""" ^ 4 E ^""^ • (2.53) 



Then, in analogy to ( 2.44| ) the thermodynamic formalism allows to estimate the r.h.s. of (2.53) for 



large Id- Together with g = 2 one thus observes that C{ld) ~ -^ log^c/i Id ~^ oo, with the effect that 
finally 

E^ {E) « A [log Ih + 1- log 2] , (2.54) 

where now Id no longer appears. However, if one wants to calculate a numerical value for T/^{E) 
the approximations employed might not be satisfactory so that the value of the constant on the r.h.s. 
of ( |2.52D cannot really be fixed by our considerations. The relation (|2.52[) then only shows that this 
constant depends on the distribution of the short periodic orbits (with lenghts < Id)- If in addition 
we introduce the explicit form Ih = '^irhd^/E for the Heisenberg length, we can obtain from ( p. 52 ) the 
following high energy asymptotics, 

T.l^{E)^\\ogE , E^oo. (2.55) 
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If we now compare (2.51) and ( p.52 ) we realize that the number variance, as a function of L, grows 
logarithmically up to approximately L ~ L^ax- Beyond this scale T,'^{L;E) ceases to increase, but 
rather oscillations caused by classical periodic orbits set in which take place about the saturation value 
S^(i?). Both the transition scale Lmax = In/ld and the saturation value depend on E: The larger 
E, the larger are both L^ax and S^(£'). Since the number variance is 'generically' well described 
by random matrix theory in the domain L <C Lmax, this is sometimes called the universality regime. 
On the other hand, the oscillations of the number variance about S^(i?) in the domain L ^ Lmax 
suggest to call this the saturation regime. This general picture has been confirmed numerically in the 



case of the zeros of the Riemann zeta function |30| and for some two dimensional classically chaotic 
systems, e.g. in [31, ^]. 

As a final remark on the number variance let us compare the results of the semiclassical analysis 
with the general observations made in section |2.1| . To this end we first have to recall the connection 
between the variables E and x: since x = N{E) ~ j^E = dE, one observes that Ih = 'iTrhd^/E ~ 



Airhydx = \fA^KA\fx. While discussing the relation ( p. 15 ) we noticed that the limit L ^ oo has to be 



considered in conjunction with x — > oo. Due to the above the latter is equivalent to the semiclassical 
limit /// — > oo. We also related the two limits in choosing L = ax'^ with some exponent < 7 < 1. 
Since Lmax = Ih/IcI = const, ^/x, the two cases L/Lmax -^ and L/Lmax -^ 00, which we were 
forced to distinguish in the semiclassical analysis, therefore correspond to < 7 < 2 and 2 < 7 < 1, 
respectively. We hence now have a means to characterize what we before alluded to as a large exponent 
7. The latter should result in a decoupling of the two random variables Nfi{y) and Nfi{y + L) on I^, 
when x ^ 00 and L = ax'^ . Now, choosing 7 larger than the critical exponent lc = \ corresponds to 
L/Lmax -^ 00 and therefore to the saturation regime of the number variance. Here S^(L; E) oscillates 
about the L-independent saturation value $]^(S). Recalling the relation ( p.l5| ), we hence conclude 
that the oscillations can only be caused by the correlation function on the r.h.s. of ( p. 15 ), whereas 



the -L-independent contribution S^(.E) comes from the first two terms on the r.h.s. Moreover, ( 2.55| ) 
implies that 

^- [^li] ~ ^ log X , X ^ 00 , (2.56) 



compare also (2.16|). So far the whole analysis concerned two dimensional systems. An immediate 



generalization to the d dimensional case reveals that Lmax ~ const. x^~d^ so that the critical exponent 
7c above which saturation takes place is in general 7c = 1 — 3- Thus, in larger dimension the saturation 
regime for the limit L ^ 00 becomes smaller. However, it could only disappear in infinite dimensional 
systems. 

The above analysis therefore leads us to conjecture that the two random variables Nfi{y) and 
^fi{y + L) on Ix become statistically independent in the limit x — > 00 when L = ax'^ with 7 > 7c. 
Once 7 < 7c, strong correlations remain that force the number variance to follow the predictions of 
the GOE. For certain two dimensional classically integrable systems Bleher and Lebowitz |^^ proved 
that indeed the saturation behaviour, with Ex{N'^j\ ~ const, ^fx^ occurs as well as the statistical 
independence for 7 > 7c = 2- '^'^^ difference in iE2;[iVrJ between classically integrable and chaotic 
systems (|2.56|) stems from the occurrence of periodic orbits as one-parameter families in the former 



case. Under certain regularity assumptions, Sarnak |36| showed that the existence of one-parameter 
families of closed geodesies in the unit tangent bundle over a surface implies Y?{L\E) » \fE for 
L ^ E ^ 00. To this end he used the trace formula in order to represent Nfi{E), and in particular 
the structure ( |1.35| ) implied by fe-parameter families which yield exponents nij = k + 1. In the example 
of the stadium billiard, where the bouncing ball orbits occur as a one-parameter family that can be 
dealt with in full detail, Sieber et al. [|^ determined the saturation value of the number variance as 
S^(i?) = const. vE, and thus confirmed the effect of non-isolated periodic orbits even in classically 
chaotic systems. 
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3 Prom Local to Global Eigenvalue Statistics 

The principal measure for the distribution of the eigenvalues of some quantum Hamiltonian that was 
discussed in section Q was the probability E{k; L) to observe k defolded eigenvalues in an interval of 
length L. As long as L is finite one thus measures correlations on a scale L, i.e., among approximately 
L eigenvalues. One is therefore dealing with local eigenvalue statistics. If one wants to include 
correlations among infinitely many levels, one is forced to consider the limit L — > oo. At several stages 
in the previous discussion we already considered asymptotics for large L, and hence the transition to 
the global scale. It is the goal of the present section to investigate this limit in some more detail. We 
will in particular emphasize the distinction as to whether one passes to the global scale within the 
universality regime L/Lmax -^ 0, or in the saturation regime L/Lmax -^ cc. It will then be noticed 
that the latter approach offers an interesting opportunity to identify fingerprints of classical chaos in 
the distribution of eigenvalues. 

3.1 The Renormalized Random Variable 

Our discussion of the local eigenvalue statistics in section ^ was based on the random variable nL{y) = 
N{y + L) — N{y) for y S /^ = [x — A.x,x + Ax]. We noticed that its expectation value reads 
JExini] ~ L, X — > oo, see ( |2.9| ). Moreover, its variance is essentially given by the number variance 
( |2.12| ) , which grows logarithmically for classically chaotic and linearly for classically integrable systems 
as L — > oo once one confines oneself to the universality regime L <C L^ax = const, ^/x, see ( |2.51| ). 
In the saturation regime L ^ Lmax it oscillates about S^. In any case, both the mean value and 
the variance of the random variable ni^y) diverge as L — > oo, x — > oo. One is therefore advised to 
renormalize the random variable so as to yield finite first and second moments in order to have a 
chance to obtain a finite limit distribution. The most simple choice 

f ^ nL{y)-L 

certainly ensures that the lowest two moments are finite, since 

1 r^+Ax ]EJnL]-L , , 



for X -^ oo, and 



IE, 



{UL - L) 



^- Wl] = W7^ = 1 • (3-3) 



In general, the ?TT,-th moments of the renormalized random variable r]L{y) can be expressed in terms 
of the moments Tj'^{L), see ( p. 11 ), of nL{y) — L, 



1 y^d) 

M^{L) := E, [C] = ^ ^x [{riL - L)""] = "^ L ■ (3.4) 

[S2(L)]T [S2(L)]^ 



Due to the normalization of i]L{y), as expressed by ( |3.2D and ( |3.3D , one is hence led to expect a limit 
distribution for L ^ oo to exist. 

The characteristic function of iji^y) can now be expressed in terms of the characteristic function 
oiriLiy), 



-^^ ^ 



J,„.(0 = Ex [e^^"-] = e" V^^ Jn„x (-/p^j ' (3.5) 
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compare (2.17). This quantity allows to characterize the distribution of the random variable r]i{y) 
through 

^VLAdr]) = TTT— / 6{rjLiy) - v) dy di] 

^^X Jx-Ax 

fX+Ax 1 r+oo 



2Ax 

r'+OO 



/ ir S'^^^y^-^'^^didydr] (3.6) 

Jx-Ax 2vr 7„oo 



1 r+oo 

Going back to ( |2.18| ), we reahze that 

y^mr)E,(vVmL)+L;L)dr^ = :^^^ j^\-^f>in^M^)^L) J^^^,[p)dpdr^ 

^ .+^^/e|(I) 
= ir , e-<^J^,,,[Odid^, (3.7) 

where we changed variables from p to ^ := pyjT!^{L) and employed (^l5|). Hence a comparison of (|3 
and (|3.7D yields 

vr,,,x{dv) ~ VWl)E, (vVWl) + L- l) dri , (3 



whenever the number variance approaches infinity, i.e., when L — > oo, x — > oo. Since we concluded 
( p.8D from the relation ( p.5| ) of the characteristic functions, the asymptotic equivalence of both sides 
of (|3.8| ) is to be understood in the sense of weak limits. This means that for any bounded continuous 
function g{r]) 

g{v) i^vLAdv) = Jim ysf(Zy / gir^) E, U^/W^) + L; l) drj . (3.9) 

As a first example, let us discuss a Poissonian random process ^^(y). Since this is stationary there 
appears no dependence of its distribution on x. The probability to find k events (i.e., eigenvalues) in 
an interval of length L is then given by 

T k 

Epoissonik; L) = — e-^ . (3.10) 

Now, it is well known and easy to see that ( 3.10| ) approaches a Gaussian for L ^ oo, 

1 (fc~L)2 1 _(k-Elnjjf_ 

since IEpoisson[nL] = L and ^'poisson(^) ~ ^- Hence, according to ( p!q ) the distribution of r]L{y) 
approaches 

1 1 2 

^VLidv) -^ ^=e"2'' dri , L ^ oo . (3.12) 



In random matrix theory one obtains a similar asymptotic result. Here the random process nL{y) 
is also stationary, after one has performed the limit of matrix dimension N —^ oo, and hence no 
x-dependence occurs either. In section 2.2 we already reported on the numerical finding (2.24). Due 
to the relation (|3.8D one therefore would conclude that ( 3.12| ) also holds for the GOE. Indeed, Costin 
and Lebowitz |39| calculated the cumulants of the distribution of r]L{y) in the GOE, GUE, and GSE 
(Gaussian symplectic ensemble) and proved that Vru^idr]) weakly converges to a Gaussian with zero 
mean and unit variance. 
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At this point we recall that the stationarity of the random process in the Poissonian as well as 
in the random matrix case makes the limit L — > cxd unique, i.e., no saturation effects and related 



phenomena occur. As we saw in section 2.3 the situation is different for spectra of individual quantum 
Hamiltonians. Now, provided one accepts that the local eigenvalue distributions of 'generic' quantum 
Hamiltonians are described by Poisson statistics in the classically integrable case and by random 
matrix theory in the chaotic case, and that this description extends to hold when passing to the global 
scale within the universality regime, the above findings suggest that 

9{v)^VLAdv) = —f^ \ 9{v)e~^'' drj , (3.13) 

-oo V^TT J-oo 

for L -^ oo, X -^ oo, L/Lfnax -^ 0, and any bounded continuous function g{r]). Thus, a central 
limit theorem would hold true for the global distribution of eigenvalues in the universality regime, 
irrespective of the type of classical dynamics. 

As can be anticipated from the discussion in section |2.3|, the situation changes drastically, if 



one passes to the global scale within the saturation regime L/Lmax -^ oo. In 2.3 we interpreted 



the saturation phenomenon in terms of the relation ( 2.15 ) for the number variance and concluded 
that all evidence seems to be in favour of an asymptotic statistical independence of the two random 
variables Nfi{y + L) and Nfi{y), y G I^, as L/Lmax — > oo. We remark that instead of choosing the 
arguments y + L and y, with y G Ix, one can obviously also consider Nfi{y) with y £ Ix and y £ Ix+L, 
respectively. In analogy to the discussion above, one should also renormalize these random variables 
and hence introduce 

Wiy) := -p^ , (3.14) 

for y £ Ix and y € Ix+L, respectively. The first and second moments of this random variable are 
clearly analogous to (|3.2| ) and (^.31) , 

]Ex[W]^0, x^oo, and lEx[W'^] = l. (3.15) 

In view of the saturation of the number variance, 

S2(L) ~ S^ ~ 2]Ex[NJi] ~ 2]Ex+L[NJi] (3.16) 

for a; — > oo, L — > oo, L/Lmax -^ oo, the random variable 



appears as the difference of the asymptotically independent random variables W{y), y G Ix, and W{y), 
y e Ix+L- Thus, instead of the limit distribution of rjL{y) as L/Lmax — > cO; one can equivalently 
consider the limit distribution of W{y) as x ^ oo. Provided the statistical independence alluded to 
above indeed holds, both limit distributions coincide. 

The quantity W{y) has a nice and simple interpretation in that it describes the -normalized- 
fluctuations of the spectral staircase N{E) about its mean behaviour N(E). Furthermore, it can 
easily be calculated from numerical data, and its distribution is readily computed. It thus appears 
that the global eigenvalue distribution in the saturation regime is most conveniently measured in terms 
of the limit distribution of W{y) for large x. 

3.2 Classically Integrable Systems 

In contrast to the case of classically chaotic systems, spectra of quantum Hamiltonians that arise 
as quantizations of integrable classical systems allow for a fairly explicit treatment. Due to the 
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semiclassical EBK-quantization scheme the behaviour of the spectral staircase N{E) can be related to 
a lattice point problem. Given an integrable classical system with d degrees of freedom, one introduces 
action-angle variables {1,0). The classical Hamiltonian function then depends only on the actions, 
H{I) = H(Ii, . . . , Id). According to the EBK-procedure, semiclassical approximations to the quantum 
energies can be obtained once one quantizes the actions as I^ = {n^ + ^)h. Here n^ S Z, or n^ S N, 
yields a quantum number, and a^ S {0, 1, 2, 3} is an appropriate Maslov index. Then 

E^FZn. ■.= H((n, + ^)h,..., (n, + ^) h) (3.18) 

defines semiclassical energies which approximate the eigenvalues of the quantum Hamiltonian H for 
small h. The EBK-spectral staircase 

N'^^^iE) := # {(m, . . . ,nrf); <^^„^ < E) (3.19) 

counts the number of lattice points (ni, . . . ^nd) inside the domain that is bordered by the {d — 1)- 
dimensional manifold defined by H{Ii, . . . ,I(i) = E. As ususal in lattice point problems, the leading 
behaviour of N^^^{E) in the limit i? ^ oo is provided by the volume of the manifold which is given 
by the condition H(Ii, . . . ,1^) < E. If one then takes the complete phase space into account by 
introducing an additional integration over the angle variables one observes that 

N''^^{E)r.^-L^I JeiE-Hil)) dldO. (3.20) 



The r.h.s. is identical to the r.h.s. of ( |1.2D because the transformation (/, 0) i— > (p, x) is canonical and 
hence preserves phase space volumes. Therefore ( 3.19D yields the actual spectral staircase at least 



to leading order as /i ^ 0. However, as discussed at length in section g the fine structure in the 
eigenvalue distribution, which is of interest in quantum chaos, is encoded in the remainder term to 
the leading asymptotics (|1.2|) of the spectral staircase. In order to analyse spectral statistics in terms 



of the lattice point problem ( 3.19| ) one therefore has to ensure that the limit distributions arising from 



N{E) and N^^^{E), respectively, coincide. For the free motion of a particle on a surface of revolution 



that obeys a certain non-degeneracy condition Bleher |55] proved this to be the case. The same fact 
was proven by Kosygin et al. [^ for the free motion on a torus with Liouville metric, which is also 
integrable. 

As a first and very simple example of a classically integrable system let us discuss a one dimensional 
billiard, i.e., the free motion on the interval [0,/] with elastic refiections from the boundary points at 
and /. In units where h = 1 = 2m the eigenvalues of the quantum Hamiltonian H = —-^ with 
Dirichlet boundary conditions read E^ = yr"^ , n € N. If we denote the integer part of a real number 
r by [r], the spectral staircase is given by N(E) = [—\/E] ~ —vE, E — > oo. As the defolded energy 
variable we thus introduce x := -\E. Then the fiuctuations of the spectral staircase read 



1 1 °° 1 

^//(^) = o - {^} = - y" - sin(27mx) , (3.21) 



vr '^ — ' n 



where {x\ = x —\x\ denotes the fractional part of x. The ^ appearing in the middle term of ( 3.2l| ) 



is chosen such that Nfi{x) fiuctuates about zero. Clearly, \Nfi{x)\ < ^, and the r.h.s. of ( 3.21| ) 
constitutes a representation of the one-periodic and piecewise continuous function Nfi{x) by its Fourier 
series. The variance of Nfi{x) can easily be determined to yield IEx[N'ii] ~ ^ for j; — > oo, so that 
W{y) = ^/V2Nfi{y). The characteristic function for the distribution of W{y) is given by 

JwAi) = 1^ r ^ e'^'^^'^ dy . (3.22) 
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In the limit x ^ oo we obtain 

Since the density p{w) of the hmit distribution of W{y) is related to the characteristic function by 
Fourier inversion, we observe from ( 3.23| ) that 




, \w\ > V-J 

Thus, W{y) is asymptotically uniformly distributed on a finite interval. The fact that p{w) = for 
\w\ > \/3 results from the bound |W^(y)| = \/T2\Nfi{y)\ < \/3; this clearly determines the support 
of the probability density p{w). The result (3.24) could also be naively anticipated from the explicit 



form and the periodicity of Nfi{x). 

A more prototypical example, however, of a classically integrable system is provided by the free 



motion of a particle on a two dimensional torus. As discussed at the end of section p.2| , the quantum 
energies for this system are exactly represented by an expression that is proportional to the positive 
definite binary quadratic form ( |2.26| ). The spectral staircase N[E) is hence exactly given by a lattice 
point problem; more precisely, N(E) is the number of points of the lattice Z,"^ inside an ellipse that is 
defined by the quadratic form Q. In this example we again choose units such that h= 1 = 2m. 

Let us now discuss the torus M.'^/C with £ = Z^ in some more detail. Here N{E) is the number of 
lattice points inside a circle of radius ^^, i.e., we are dealing with the circle problem. Employing the 
Poisson summation formula, one can readily demonstrate that for E ^ oo 

E Af:^ cos (kV^^ - ^) +0(e--A . (3.25) 




2vr3 „g^ mi f,^^ k2 \ 4 



Since ^^ ^ {nf + 722)" * diverges, one cannot estimate the sum on the r.h.s. of ( 3.25| ) in a simple way 
as a function of ii^; to obtain improved bounds for this sum still constitutes an active area of research 
in analytic number theory. A novel perspective on the circle problem, which inspired much of the 
work on the global distribution of quantum energies, was introduced by Heath-Brown |^]. Instead of 
estimates of the remainder term in the circle problem, he considered its distribution. As we will soon 



realize, this corresponds to the limit distribution of W{y). According to ( |3.25| ), where the first term 
on the r.h.s. defines N{E), the defolded spectral variable reads x = ^. Thus 

Nfiix) = (^-y <^{x) + O (^x~^^ , x^cx), (3.26) 

with 

m square free 

where 

fm{t) := -\ Y. -^ ^°^ f ^vrfcf - ^\ . (3.28) 

vrm4 ^^-^ k^ \ * / 

It is furthermore known, see for example [^, that iE'^[A^?J ~ const, ^fx., x -^ 00, so that W{y) = 
const, ^{y). Apart from the normalization of the variance, the limit distributions of W{y) and $(y) 
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therefore coincide. Notice that the representation (3.25) of the spectral staircase can be interpreted 



in terms of a semiclassical trace formula. To this end one observes that \Jvi\ + n^, with (ni,n2) G 
Z'^yKO, 0)}, yields all lengths of periodic orbits on the torus. In the same spirit, the sum in ( [3.27] ) can 
be identified as extending over all lengths \pm of primitive periodic orbits; their repetitions are then 
accounted for by the sum over k in ( ^.28 ). 

The functions fmit) in ( |3.28| ) are clearly continuous and periodic with period 1. For such functions 
Heath-Brown proved |5^] that whenever 71,..., 7^ are real numbers linearly independent over Q, 
then fmii'fit), ■ ■ ■ , frnkilkt), with < t < T, become statistically independent in the limit T -^ 
00. If one now chooses 7^ = \/m, m square free, any set {7mi 7^ • • • 7^ 7mfe} is indeed linearly 
independent implying that the random variables fm{t\/rn), < t < T, as they appear in ( 3.27] ), 
asymptotically become independent for T — > 00. A further analysis reveals P8| , p8|| that under fairly 
general assumptions about a set of continuous real valued periodic functions 61 (i), b2{t), . . . with 



Jo 



bn{t) dt = 



00 „i 
and 2. / bn{t) dt < co 
„=i -^0 



(3.29) 



the series F{t) = ^^ bni'Jnt), < t <T, has a distribution that weakly converges to the distribution of 
the random series YlnbriiOn)-, where the 0„'s are independent random variables uniformly distributed 
on the interval [0, 1]. The latter fact arises because the periodicity of the functions hn{t) results in 
7„,t to be taken mod 1. As t varies over the interval [0,T], the values of 7„i mod 1 become uniformly 
distributed on [0, 1] as T — > 00. 

We remark that for the following results to hold the second condition in ( 3.2S| ) becomes essential; it 
roughly means that the maximum of hn{t) must decrease sufficiently fast as n ^ 00. Then a sufficient 



condition for the existence of the limit distribution is for example given by |57, 28 1 



lim lim sup — 

iV->oo T^oo -^ 



min < 1, 



F{t) - Y. ^"(^"*) 



n<N 



> dt = 



(3.30) 



the infinite sum 
it was moreover 



This means that F{t) can be represented, in the sense indicated in ( |3.3[1| ), by 

Yln^riiln't)- Under further assumptions on the behaviour of the functions 6„(i) 

proven ||2^ that the limit distribution has a density p{x) which behaves asymptotically as e"*^'^' for 

X — > ibcxD. For the circle problem it was verified |p8| that the functions fmif) nieet all requirements, 

and the exponent was estimated from above and below hy p± = 4ibe for any e > 0, respectively; thus 

the limit distribution is not Gaussian. 

We are interested in the limit distribution of W{y) for the eigenvalues of the quantum Hamiltonian 
on the torus M^/Z^, and hence in the limit distribution of $(y), see ( p. 26 ). By the change of variables 

,) and then 



7ri we can relate this problem to the distribution of the random series YL^ fn 



apply the results of [^, 28 1. Therefore 

1 



lim 

x^oo 2Ax 



-Ax 



9iW{y))p 



-Ax 



2Ax 



X 1\ , 
+ -^]dy 



+00 



g{w)p{w) dw 



(3.31) 



for any piecewise continuous bounded function g{w) and any probability density p{t), t E [0, 1]. Here 
p{w) is the non-Gaussian probability density mentioned above. The freedom to choose a weight p(r) 
allows to change variables from, say, y to t and thus to compensate for the Jacobian then arising. 

The existence of a limit distribution ( p. 31 ) has by now been established for a large class of two 
dimensional classically integrable systems, see |5^] for a review. This includes the free motion on 
more general fiat tori, with and without Aharonov-Bohm type magnetic fluxes p8| , |60| ], on surfaces of 
revolution |5^], and on tori with Liouville metrics |56, 61]. All these cases lead to limit distributions 
with densities p{w) that approximately decay as e~'^'^' for w — > ±00. 

In a further class of examples for two dimensional classically integrable systems the limit distri- 
butions of W{y) were explicitly determined by Schubert [^2[. He considered the geodesic motion on 
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Zoll surfaces, i.e., surfaces with the topology of a sphere which are endowed with a metric all of whose 
geodesies are closed, and all lengths of periodic orbits are multiples of a fundamental length. The 
simplest example for a Zoll surface is the sphere with its usual round metric; however, non-trivial 
deformations thereof exist. Schubert showed that Nfi{x) = ^/xQ{^/x), where Q{t) is an almost pe- 



riodic function of Besicovitch-class B^, see [55] for explanations. The trace formula then allows to 
represent @{t) by an L^-convergent Fourier series that is, apart from a possibly non-trivial Maslov 



phase, identical to ( 3.21 ). Therefore, W{y) = Q{^/y) for Zoll surfaces has the same limit distribution 
(3.24) as a one dimensional billiard. This seemingly surprising coincidence is caused by the fact that 
in both cases the classical side of the trace formula is given by a sum over multiples of one fundamental 
length. The extra dimension of a Zoll surface thus is only reflected in the overall growth of Nfi{x) by 
^/x. Notice also the difference to the case of tori where the overall growth of Nfi(x) is given by x^ . 

As a final remark we observe that the decay of the propability density piw) as w ^ ±00 is 
connected to the type of convergence of the periodic orbit sum representing Nfi{x). For the one 
dimensional billiard as well as for Zoll surfaces the series converges in L^-norm, which is considerably 
strong. In contrast, one obtains only the weaker form of convergence ( ^.SOj) for flat tori. We hence 
conclude that the weaker the convergence is, the weaker the decay oi p{w). This observation will be 
important for classically chaotic systems. We also notice that in all examples of classically integrable 
systems studied so far the distribution of W{y) has a limit as a; — > 00 that is non-Gaussian; thus in 
the integrable case no central limit theorem for the global distribution of eigenvalues in the saturation 
regime is valid. 

3.3 Classically Chaotic Systems 

In analogy to the situation for classically integrable systems that has been described above, the global 
distribution of eigenvalues of quantum Hamiltonians with chaotic classical limits will now be studied by 
means of the limit distributions of W{y), y G Ix, as x — > cxd. Above we realized that these distributions 
can be obtained from periodic orbit sums for Nfi{x). In all of the integrable cases discussed we noticed 
that Nfi{x) = x"' (p{x), where a is some appropriate power that describes the overall growth of the 
fluctuations Nfi{x) of the spectral staircase as x — > 00, and (j){x) is a sum over the (primitive) lengths 
of classical periodic orbits and their multiples; one in particular finds that lEx{N'^i\ ~ const. x'^"', 
X — > 00. The periodic orbit sum for (^(x) was found to share some controllable convergence properties; 
it either was a Fourier series or, more typically, an almost periodic function whose representation by a 
trigonometric series converges in some Besicovitch norm. Now focussing on the chaotic case, we realize 
that a representation of Nfi{x) in terms of a sum over (primitive) periodic orbits, see for example 
( |1.63| ), results in a considerably worse convergence. Due to the exponential proliferation ( 1.64 ) of the 



number of periodic orbits, a Gaussian smoothing was required in (1.63) in order to achieve a finite 
sum. Moreover, no power of x or logx drops out explicitly in front of the periodic orbit term that 
could directly indicate the overall increase of Nfi{x) as x — > 00. However, due to ( p.56|) lEx[N'ji] is not 
bounded but rather grows logarithmically. Hence the overall increase of Nfi(x) is implicitly contained 
in the periodic orbit term itself and has to be uncovered by careful estimates of the sum. This remark 
should indicate the non-trivial problems that go along with a quantitative analysis of the distribution 
of W{y) for classically chaotic systems. 

A strong hint at what has to be expected in the chaotic case is provided by a comparison with the 
distribution of W{y) for the non-trivial zeros of the Riemann zeta function. At the end of section ^^ 
we pointed at the close analogy between the Selberg zeta function with its spectral zeros on the one 
hand, and the Riemann zeta function with its non-trivial zeros on the other hand. The asymptotics 
( |1.86| ) of the counting function for the Riemann zeros now suggests to introduce the defolded variable 
^ log ^. A 'semiclassical' analysis in the spirit of section p.3| then reveals that Ih = const, logx 



X 



has to be chosen so that also Lmax = const, logx. According to ( p. 54 ) one therefore concludes that 



lEx[N'ji] ~ 2:^^ log log 2;, X — > 00. In order to arrive at this result one has to take into account 
that lengths of periodic orbits correspond to logarithms of primes. Since the latter obviously have 
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multiplicity one, the semiclassical analysis requires to choose g = 1, as in the case of classically chaotic 
systems without time-reversal invariance. The latter point is in accordance with the conjectured -and 
numerically confirmed- GUE behaviour of the local statistics of the Riemann zeros J6^. Already in 
1946 Selberg proved 



lE^ 



N: 



2k 



fl 



SO that the even moments of 



Ml 



W{y) 



-^ log log X 



Nfiiy) 



-^ log log X 



+ 



(log log x)'^ 2 



y e 4 



(3.32) 



(3.33) 



converge to the even moments of a standard Gaussian as x ^ oo. It was later noticed that an 
application of Selberg's method to the odd moments of Wiy) shows that these asymptotically vanish, 
see for example [|65|| . The distribution of W{y)^ y S /x, therefore weakly converges to a Gaussian of 
zero mean and unit variance as x ^ oo. Selberg |6^ extended this result to the corresponding limit 
distributions for a whole class of number theoretical zeta functions with functional equations. Since 
the zeta functions covered by Selberg's theorem are required to possess representations as ordinary 
Dirichlet series, the case of the Selberg zeta function 2'(s) is not included. But based on the analogy 
between the zeros of the Riemann zeta function and quantum energies of a classically chaotic system 
as well as on numerical evidence, in p^ the following conjecture was established: If the classical 
dynamics are chaotic, with only isolated and unstable periodic orbits, and a scaling behaviour as in 
( 1.401 ) holds, the limit distribution of W[y) is a standard Gaussian. In analogy to (3.31) it hence is 
conjectured that 



lim — - — 



x+Ax 



9{W{y))p 



z— Ax 



y — X 1 

2Ax "^2 



dy 
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+00 



g{w) e 



dw 



(3.34) 



for any bounded continuous function g{w) and any probability density />(r) on [0, 1]. The restriction 
of the conjecture to chaotic systems with scaling ( [1.40 ) stems from the analogy to the distribution 
of the zeros of zeta functions. Exactly in the scaling case one can directly encode the periodic orbit 
contribution to the spectral staircase in the complex phase of a semiclassical zeta function which is a 
function of the scaling variable E'^, compare (1.40). Nevertheless, we expect the conjecture to hold 
also in more general situations, with the essential requirement that all periodic orbits be isolated and 
unstable. We also remark that the conjecture includes the case of arithmetic quantum chaos, where 
the local distribution of eigenvalues strongly deviates from the 'generic' situation for classically chaotic 
systems in that it rather shows a Poissonian behaviour than a GOE one ||33| , |34| . Concerning the global 
distribution of eigenvalues the differences between arithmetic and 'generic' systems only show up in 



the normalization of W{y) because in the arithmetic case one obtains ]Ex[N% 



\fx 



log a; ' 



OO, 



see [29|. A preliminary version of the conjecture ( p.34 ), which concerns W{\j) evaluated only at the 
eigenvalues, is contained in ||68|| . 

In ||6^ the distribution of A^// (y) has been computed numerically for two chaotic billiards in the unit 
disc ( |1.67| ) with the metric of constant negative curvature, one of them showing arithmetic quantum 
chaos, as well as for one integrable billiard. The normalization leading to the random variable VK(y) 
was not taken into account. But since in numerical calculations the spectral interval Ix from which 
y is taken has to be kept finite, the non-existence of a limit as x — > oo does not become apparent in 
numerical data. The numerical results of 1 6£ ] show distributions that are somewhat closer to Gaussians 
in the two classically chaotic examples than in the classically integrable case. A detailed numerical 
analysis of the distribution of W{\j) is presented in |7^]. There six chaotic systems are compared with 
one pseudo-integrable and five integrable systems, and a clear evidence for Gaussian limit distributions 
to arise only in the classically chaotic cases is obtained. 
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The influence of periodic orbits on the hmit distribution of W{y) entered our above discussion only 
indirectly through the representation of Nfi{x) by the complex phase of the Selberg zeta function or 
its semiclassical analogue. If one desires to uncover the role played by unstable periodic orbits in a 
more direct fashion one has to express Nfi{x) in terms of the Gutzwiller -or Selberg- trace formula. 
For this purpose let us again concentrate on the situation covered by the trace formula ( |1.48| ), or 
more specifically by the Selberg trace formula ( |1.7lD . Then Nfi(p), where p denotes the appropriate 
momentum variable, can be represented by N^jiip) for e — > 0, see ( 1.63| ). Alike in the semiclassical 
analysis of section 2.3 we now rewrite the sum over primitive periodic orbits 7p as a sum over their 
respective lengths li < I2 < h < ■ ■ ■, thereby introducing the multiplicities Qn of /„, n = 1,2,3, .. .. 
That way Nfi{p) turns out to be a sum over the functions 



""W^=tEi: 



J r^fJ'-n 



k=l 



det(M^ 



sm(2iTkt) 



(3.35) 



evaluated at t 
zero mean, 



Moreover, 



cr^ 






From ( 3.35 ) we obtain that a„(t) is a one-periodic continuous function with 



a„(t) dt = 



al 



00 



1 



i«4t)rdi = 7S>.T^ 
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00 



(3.36) 



(3.37) 



If we now employ the periodic functions a„(i) in order to express Nfi{p), we have to consider gf^^n 
mod 1. In case the lengths In of primitive periodic orbits are linearly independent over Q, 2^1^" mod 
1 becomes uniformly distributed on the interval [0, 1] as p ^ 00 or /i ^ 0. Furthermore, due to the 
result of 1 57] already alluded to above for integrable systems, the linear independence of /„^, . . . ,ln^ 
implies 



lim — - — 

p^oo 2Ap 



P+Ap 



'"rii 



p—Ap 



P 
2-nfi 



"1 



'•rik 



P 

27rh 



Irii 



— = TT 



2TTh 



,{t)dt , 



(3.38) 



where the interval \p — Ap, p + Ap] corresponds to the spectral interval Ix upon changing variables 
from the defolded energy variable y to the momentum variable p' . Thus Nfi(p') appears as a sum of 
the asymptotically independent random variables ani-^^ln)- 

When comparing the present situation with the classically integrable case one first of all notices 
that the required linear independence of the primitive lengths is not known to hold in a single classically 
chaotic example. Secondly, and more importantly, the sum 



00 00 ^]^ 

n=l n=l"^0 



dt 



(3.39) 



00 



can be shown to diverge. To this end we remark that obviously the leading behaviour of a"^ as n 

is determined by the contribution from /c = 1 in (3.37). If we furthermore recall the discussion of the 

r.h.s. of ( p. 53 ), we obtain that in case g = 2 

1 ^^ o „ Q ^^ .. 1 



27r2 



Y.9n' 
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27r2 
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log I 
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CO 



(3.40) 



Thus, if we consider the random variables an{On), where the ^^'s are independent random variables 
uniformly distributed on the interval [0, 1], their variances a^ sum up to 

N 



-Bat := V] o-^ ~ ^ log In '■ o l°g 



n=l 



- log(2iVlog27V) 
r 



-2 log log iV , 



N 



CXD . 



vr 



(3.41) 
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The asymptotic relation I 



jy ~ 1 log(2A^log2A^), N —>■ CO, leading to ( |3.41| ) immediately derives 
from ( 1.64 ). For the arithmetic systems with their exceptionally large multiplicities of lengths [^] 
Bn diverges more strongly than in the 'generic' case covered by ( 3.41] ). In the classically integrable 
situation we noticed that both conditions (3.29) ensured the existence of a limit distribution and 
allowed for conclusion on its properties. In the present situation the second condition is violated, 



which indicates bad convergence properties of the sum ^^ an{t). However, it is suggested in |59] that 
due to the divergence of ( ^.39 ) one should employ the Lindeb erg- Feller version of the central limit 
theorem, see for example |7l| . This states that the distribution of the random variable 



iN-.-- 



E 



N 
n=l 



(3.42) 



converges as A^ ^ oo to a Gaussian with zero mean and unit variance. 

We now again adopt the point of view, already employed in section ^, that the periodic orbit 
sum for Nfi{p) is dominated by the terms with lengths In ^^H = 4:'irphd. If we then denote the index 
corresponding to Ih by Nh , Inh = Ih, ^^ first of all notice that Nfi (p) is essentially given by 



Nh 



P 

2TTh 



l 



(3.43) 



We now also cut off the sum of the variances o"^ that refer to the individual terms in ( 3.43| ) at Nh, 

and hence obtain 

Nh 



1 



1, 



Bnu = E ^" ~ ^ l°g ^^ ~ 2 ^^^^^ ~ ^" ^'^ 



(3.44) 



n=l 



see ( 3.41| ), (|2.54 ), and ( 2.56| ). The relation ( p. 44 ) is remarkable since it states that in the semiclassical 
limit the true variance of Nfi{p), which appears on the very r.h.s., is asymptotic to the sum of the 
variances o"^ when this is cut off at Nh. Therefore the two assumptions: statistical independence of 
the a„'s and truncation of the periodic orbit sum at Ih, merge into a coherent picture. We conclude 
that thus ^Nh essentially yields W{y), 



Wiy) 



Nfijy) 



-Nh 



Cnh 



V^Nh 



(3.45) 



and that the limit distributions coincide. Here y £ Ix, and x corresponds to Ih via Ih = 'inhVdx. 
Notice that the semiclassical limit //^ — > oo has to be performed in (3.45) as x ^ oo on the l.h.s., and 
as Nh ^^ oo on the r.h.s. We therefore expect the limit distribution of ^at as A^ ^ oo, see ( 3.42| ), to 
reproduce the limit distribution of W{y), y £ Ix, as x ^ oo. At this point we want to stress that the 



above considerations do not constitute a proof of ( 3.34| ), since the assumptions made at various stages 
remain to be verified. However, we believe that all evidence is in favour of the conjecture. 

The above discussion of the limit distribution of W{y) in the classically chaotic case was restricted 
to systems with only isolated and unstable periodic orbits. Only for those the zeta function analogy 
works and the periodic orbit expression for Nfi{x) has the form required by ( p. 35 ). In order to illustrate 
the necessity of the condition on the periodic orbits let us discuss the stadium billiard in some detail. 



The classical dynamics are known to show the K-property |72|, and hence the stadium billiard qualifies 
by all means as a chaotic system in the usual sense. However, not all periodic orbits are isolated and 
unstable because of the presence of the bouncing ball orbits that form a one-parameter family of 
periodic orbits. In addition to the contribution of the isolated and unstable orbits the Gutzwiller 
trace formula for the stadium billiard therefore contains a different term caused by the bouncing ball 
orbits. Thus the spectral staircase expressed in the momentum variable reads 



N{p) = N{p) + NfiM + Nbbip) , 



(3.46) 
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where Nfi^u{p) is the ususal contribution of the unstable periodic orbits. In |54] an expUcit expression 
for the contribution of the bouncing ball orbits is derived, 

iVfe6(p) = ^^=P^ V — COS ( 2anp - ^ ) + O L-^\ . (3.47) 

2V7r-^a ^ n2 V 4 y \ / 

Here h denotes the length of the two parallel edges of the billiard domain, and a is the radius of the two 
half-circles. In fact, ( p. 471 ) corresponds to the desymmetrized stadium billiard whose domain consists 
of a quarter of the stadium domain, see |p3 for details. In terms of the defolded energy variable 
X = di? the variances of the two periodic orbit contributions in (|3.46D read 

^^ [^A J ~ A log ^ ^^d ^- [^'fc] ~ ^ f^^l ^^ , X ^ OO . (3.48) 



The first term is just (p. 56), whereas the second term was determined in |54|. Hence 



lE^ [(iV/,,„ + iVfeb)'] ~ -^T^ ^^ ' ^ ^ OO > (3-49) 

so that W{x) is asymptotically completely determined by the bouncing ball contribution 



, , N^h{x) / 2 ^ 1 / /J 37r\ 



^A^L\ 



,21 VC(3)f-ni V Vd 4 7 



n=l 



Thus VFfeft(x), considered as a function of -^x, is asymptotic to a continuous periodic function with 
Fourier series given by ( p. 50 ). Since X]n^~^ < oo we can moreover estimate 




|T^66(x)|<. -— V— = W— -CU ■ (3-51) 



This finite bound implies that the density p(w) of the limit distribution of M^(y) for the complete 
desymmetrized stadium billiard vanishes outside the finite interval [— c, +c], where c denotes the con- 
stant on the r.h.s. of ( ^.51 ). The limit distribution is therefore clearly not Gaussian, and this violation 



of ( |3.34| ) is solely caused by the bouncing ball orbits. We remark that this observation is not restricted 
to the stadium billiard, but will occur in all systems with families of periodic orbits. These always 
yield variances of their contributions to the spectral staircase that exceed the logarithmic contribu- 
tion ( p. 561 ) of the unstable orbits. By dividing through the rms fluctuations of the spectral staircase 
in the definition of W{y\ the influence of the otherwise dominating unstable orbits hence vanishes 
asymptotically. 

3.4 Summary of the Global Statistics 

After having discussed various aspects of the global distribution of eigenvalues for classically integrable 
as well as chaotic systems in some detail we now want to summarize the emerging perception of the 
global statistics. But before we again stress some points that we already remarked before at several 
places. 

1. Statements about local as well as global eigenvalue statistics always refer to 'generic' systems in 
a given class. A precise definition of what this means has to be given in each particular case. 
A typical example would be some statement concerning eigenvalue statistics on tori, see ( 2.26| ), 



where one would have to specify the set of parameters A^B^C to which the statement should 
apply. It seems to be almost impossible, however, to give a general characterization of a 'generic' 
classically integrable (or chaotic) system in a unified manner. 
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2. Concerning the local statistics of eigenvalues it is assumed that it can 'generically' be described by 
the appropriate random matrix ensembles for classically chaotic systems ||4l|] , and by Poissonian 
statistics for classically integrable systems 



3. We furthermore assume that the description via random matrix ensembles or Poisson statistics 
extends to hold in the limit L — > oo once this is performed in the universality regime L/Lmax — *■ 0. 
We base this assumption on the following observation: The local statistics of eigenvalues emerges 
as the limit of x ^ oo of, say, Ex{k; L) for finite L. Passing to L ^ oo within the universality 
regime appears to leave the analogy to, say, the GOE in the number variance unchanged, see 



In the course of our discussion of eigenvalue distributions it became clear that the random variable 
VLiy), y ^ Ix, appears to be most suited for a unified treatment of spectral statistics on all scales. After 
the limit x — > oo has been performed, the distribution of rjL with finite L yields the local statistics. 
According to (|3.8| ) and the assumption stated in 2., then 

i^rjAdT]) = V^^iL) E (vV^HL) + L;L) dr] , (3.52) 



where E{k; L) is either given by random matrix theory or by Poisson statistics. 

The transition to the global scale now proceeds under the assumption in 3. We hence conclude 
that the distribution of rii{y), y £ I^-, converges for x -^ oo, L — > oo, L/Lmax ^ to the same limit 
distribution as the respective distributions of r]L in the random matrix ensembles or for a Poisson 
random process when L — > oo. Since the latter limit distributions are standard Gaussians, see [^] 
and ( p. 11 ), we therefore expect the central limit theorem ( 3.13| ) to hold for all 'generic' systems, 



independent of their respective classical limits. It thus appears that the global spectral statistics in 
the universality regime are not sensitive enough a measure of eigenvalue correlations in order to yield 
information about specific features of a quantum system. 

Finally we have to consider the limit x -^ oo, L -^ oo, L/Lmax -^ oo, i.e., the transition to the 



global scale in the saturation regime. In section 3T we argued that the limit distributions of rji^y) 



y £ Ix, then coincide with the limit distributions of W{y), y E I^;, as x — > oo. This conclusion 
was drawn from the apparent asymptotic independence of the random variables Nfi(y), y G Ix, and 
^fi{y)^ y £ Ix+L, as L/Lmax -^ c>o, which was discussed at the end of section ^. The principal 
observation then made in sections W^ and B^ was that the limit distributions of W{y) allow to 
distinguish classically integrable systems from classically chaotic ones with only isolated and unstable 
periodic orbits. Only the latter yield Gaussian limit distributions, whereas in all other cases the 
distributions of W{y) are found to converge to a non-Gaussian limit. 

We are now in a position to summarize the above findings in the following generalization of the 



conjecture ( 3.34 ) introduced in [67|. Considering the limit distribution of the random variable rji^y), 
y E Ix, we claim that 

f + OO /■ + 00 

lim / g{r])vr,L^x[dr])= / g{ri) P{ri) dr] , (3.53) 



X — >oo 

L-,oc ■J —CO J —CO 



for any bounded continuous function g{rj). The limit distribution P{rj)dr] shall have a density P{ri) 
that depends on the way the double limit is performed, as well as on the type of the classical dynamics. 

1. L/Lmax -^ 0: For 'generic' systems, be they classically integrable or chaotic, one obtains 

P(r]) = ^e-h\ (3.54) 



2. L/Lmax — *■ oo: The limit density P{ri) depends on the type of the corresponding classical system. 
When the latter is integrable (or chaotic with families of periodic orbits) P{r]) is non-Gaussian. 
Most probably it will decay faster than a Gaussian, 

lim P{rj) e^''" = . (3.55) 

r;— >itoo 
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In certain cases it can even have a compact support. If the classical dynamics are chaotic, 
with only isolated and unstable periodic orbits, then a standard Gaussian as in ( |3.54 ) is to be 
observed. 

In the field of quantum chaos so far mainly the local distribution of eigenvalues was investigated; 
in particular the distribution of nearest neighbour level spacings played an eminent role. The results 
as described by the random matrix/Poissonian conjecture, however, seem to be somewhat counter- 
intuitive. A Poissonian random process, which describes classically integrable, i.e. regular, systems, is 
characterized by the absence of any correlations. On the other hand, the behaviour of classically chaotic 
systems shows strong correlations. In this respect the situation concerning the global distribution of 
eigenvalues in the saturation regime seems to be more confirm to intuition. The Gaussian limit 
distribution emerging for W{y) in the chaotic case is the most random of all possible probability 
distributions that have a density and are of a fixed common variance. If we introduce the spectral 
entropy 

r+ca 

£[p] := — / p{w) logp{w) dw , (3.56) 



which provides a quantitative measure of a mean unlikelihood for W{y) to have a certain value, the 
distribution with maximal entropy characterizes the most random spectral fluctuations. It is well 
known that under the constraint of a fixed variance, £[p] is maximized by a Gaussian of zero mean. 
Therefore, quantum spectra that have a Gaussian limit distribution for W{y) or ijiiy) possess, on a 
global scale, the most random fluctuations. This seems to be a convenient characterization of quantum 
chaos. 
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